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Preface 


Preface 


“Essential Electromagnetism” and “Essential Electrodynamics” (also to be published by Ven- 
tus) are intended to be resources for students taking electromagnetism courses while pursuing 
undergraduate studies in physics and engineering. Due to limited space available in this series, 
it is not possible to go into the material in great depth, so I have attempted to encapsulate 
what I consider to be the essentials. This book does not aim to replace existing textbooks 
on these topics of which there are many excellent examples, several of which are listed in the 
bibliography. Nevertheless, if appropriately supplemented, this book and “Essential Electro¬ 
dynamics” could together serve as a textbook for 2nd and 3rd year electromagnetism courses 
at Australian and British universities, or for junior/senior level electromagnetism courses at 
American universities/colleges. 

The book assumes a working knowledge of partial differential equations, vectors and vector 
calculus as would normally be acquired in mathematics courses taken by physics and engineering 
students. However, very brief introductions to vectors, vector calculus and index notation are 
given in the appendices. Some of the mathematical derivations have been relegated to the 
appendices, and some of those are carried out using index notation which is briefly introduced 
in Appendix D, but elsewhere in the book manipulation of equations involving vector differential 
calculus is done using standard vector calculus identities (also given in the appendices). 

The presentation of the subject material in this book is conventional, starting with forces 
between charges, electric held, Coulomb’s law, electric flux, Gauss’ law and the electrostatic 
potential. Chapter 2 is on Poisson’s and Laplace’s equations, and their solution. The method 
of images is applied to simple examples in plane, cylindrical and spherical geometry. Laplace’s 
equation is solved analytically in Cartesian coordinates for the cases where the boundaries 
are orthogonal planes, and in spherical coordinates where the boundary surface is a sphere; 
these being the most commonly-encountered problems involving Laplace’s equation. Solution of 
Laplace’s equation in cylindrical coordinates is not included, but may be found in most applied 
mathematics texts. A numerical finite-difference method is described for solving Laplace’s 
equation. This method is applicable whether or not the boundary conditions allow analytic 
solutions. 

Chapter 3 is on the multipole expansion of the electrostatic potential, up to and including the 
quadrupole term. Chapter 4 is on macroscopic and microscopic dielectric theory, starting with 
the polarisation held, polarisation charges, Gauss’ law and the displacement held, susceptibility 
and permittivity, and the boundary conditions on the electric and displacement helds. Orien- 
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tational polarisability of molecules, electronic polarisability of non-polar molecules and ionic 
polarisability of crystals are discussed. Finally, the Clausius-Mossotti formula which connects 
the microscopic and macroscopic theories is derived. 

Chapter 5 deals with the magnetic field, its causes, magnetic forces and the magnetic flux. The 
magnetic vector potential is derived from the Biot-Savart law, and a derivation of Ampere’s 
law follows, together with examples of its use in solving problems in magnetostatics. Finally 
the multipole expansion of the vector potential up to the dipole term is performed, and mag¬ 
netic dipoles are discussed together with their magnetic field and the torques and forces they 
experience in a magnetic field. 

Chapter 6 focuses on the magnetism of materials, starting with the magnetisation field, mag¬ 
netisation currents and their inclusion in Ampere’s law, and the introduction of the auxiliary 
field H. Next, susceptibility and permeability are defined, and the boundary conditions on the 
magnetic and auxiliary fields are derived. Finally, the causes of magnetisation are discussed, 
namely orientational polarisability of atomic magnetic dipoles associated with unpaired electron 
spins causing paramagnetism, and the quantum mechanical exchange interaction in some ma¬ 
terials giving rise to unpaired electron spins lining up in the same direction as in neighbouring 
atoms, and leading to ferromagnetism. 
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Each chapter is followed by several exercise problems, and solutions to these problems are to 
be published separately by Ventus as “Essential Electromagnetism - Solutions”. I suggest you 
attempt these exercises before looking at the solutions. 

I hope you find this book useful. If you find typos or errors I would appreciate you letting me 
know so that I can fix them in the next edition. Suggestions for improvement are also welcome 
- please email them to me at protheroe.essentialphysics@gmail.com. 

This book was mainly written in the evenings and I would like to thank my family for their 
support and forbearance. It is dedicated to the memory of my parents, who stimulated my 
interest in science. 

Raymond John Protheroe, 

School of Chemistry and Physics, The University of Adelaide, Australia 
December 2012 
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1 Electrostatics 

Learning objectives 

— To learn the concepts of electric charge, charge density and the electric field, and how 
they are related through Coulomb’s law and the principle of superposition as applied to 
forces. 

— To understand the concepts of electric field lines, and electric flux as the surface integral 
of the electric field. 

— To be able to derive Gauss’ law in integral and differential form, and be able to use it to 
obtain the electric field from symmetrical charge distributions. 

— To know why some materials are conductors and others are insulators (dielectrics), why in 
electrostatics the electric field is zero inside a conductor, that any charge present resides 
on the surface of a conductor, and that the electric field just outside a conductor is normal 
to the surface. 

— To understand what the electrostatic potential is, how it may be obtained from the electric 
field, and that the electric field is the negative gradient of the potential. 

— To be able to find the work done to assemble a charge distribution, to understand that 
the potential energy is stored in the electric field, and to be able to find its energy density. 


1.1 Electric charge, field and flux 


In 1785 French physicist Charles-Augustin de Coulomb (1736-1806) published his law, which 
may be stated as: the force exerted on point charge Q by point charge q a distance R away 
has magnitude proportional to Qq/R 2 and is directed along the line between the two charges. 
American statesman and scientist Benjamin Franklin (1706-1790) was responsible for naming 
the charge we now associate with electrons “negative” and the charge we now associate with 
protons “positive”, and was the first to state the law of conservation of charge. Electric charge 
is measured in coulombs (C). 

We may consider the force on charge Q as being due to the electric field E of charge q acting 
on charge Q, he. F = Q E. The electric field is defined as the force per unit test charge (unit: 
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N C 1 or Vm x ); a test charge being an infinitesimal point charge placed where the electric 
field is to be evaluated. From Coulomb’s law, the electric field of a point charge is 


E(r) 


qH 

4t re 0 i? 2 


( 1 . 1 ) 


where R = (r — r'), r' is the source point and r is the field point as in Fig. 1.1(a), and 
£o = 8.85418782 x 10 -12 C 2 N -1 m -2 is the permittivity of free space; more usually its units 
are given as (F m -1 ) where F stands for the farad which is the unit of capacitance, and is 
named after English chemist and physicist Michael Faraday (1791-1867). 

(a) (b) 




Figure 1.1: Geometry for discussing (a) Coulomb’s law for point charge q , and (b) Coulomb’s 
law for charge density p( r). 


Since the electric field is the force per unit test charge, it obeys the principle of superposition 
just like forces do. Hence, for a continuous charge distribution with volume charge density p( r) 
(C m -3 ) the electric field is (Fig. 1.1b) 


= <“> 

We may think of the electric field of, for example, an isolated positive point charge (Fig. 1.2a), 
in terms of field lines which emerge isotropically from positive charges (Fig. 1.2b). The number 
of field lines per unit area crossing a concentric sphere of radius r is then nj (47rr 2 ) oc 1/r 2 where 
n oc q is the number of field lines originating at q. At any location, the electric field’s direction 
is tangential to the field line at that location. 
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Figure 1.2: (a) Electric field vectors due to a positive charge, where the length of an arrow 
is proportional to E at the location of its tail, (b) Electric field lines starting at an isolated 
positive charge, (c) Electric field lines ending on an isolated negative charge, (d) Electric field 
lines starting on positive charge +q and ending on a nearby negative charge —q. (e) Electric 
field lines starting on two nearby positive charges +q repel each other, (f) Field lines for the 
case of two opposite charges of unequal magnitude. "S" represents a closed surface surrounding 
all charges present. Note that in these, and subsequent, figures the point charges are located 
at the centres of the white discs on which their charge is written. (The field line tracing for all 
the images in this book which show electric or magnetic field field lines has been done using 
a numerical scheme based on the 2nd order Runge-Kutta method for initial value problems 
described in Appendix F.) 


1.1.1 Electric flux 


The electric flux through a surface S is defined by 

$ E = [ E-dS = f E n dS (1.3) 

Js Js 

and this is proportional to the net number of field lines passing out through S', i.e. the number 
going out minus the number going in (Figs. 1.2b-c). In the flux integral the vector surface 
element dS is the the scalar surface element dS multiplied by the normal unit vector n at the 
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surface, i.e. dS = n dS. In the case of a closed surface the direction of dS is always defined to 
be outward from the volume, but otherwise there is an ambiguity and the direction must be 
specified. 

For a spherical surface of radius r centred on the origin as illustrated in Fig. 1.3(a) 

dS = (r dQ){r sin 0 d(j>) — r 2 sin 6 d6 d<j) — r 2 |d(cos0)| d<j) (1.4) 

and n = r. Thus, for charge at the origin, the flux outward through the sphere of radius r 
centred on the origin is 

i E ' iS = L = T (L5) 

Of course, we could obtain this result more easily by exploiting the spherical symmetry of the 
problem — the surface integral is simply the product of E — q/(4:7T£or 2 ) with the total area of 
the sphere S — \n rr 2 . 
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Figure 1.3: (a) Surface element dS of a sphere of radius r centred on the origin, illustrating 
spherical coordinates (r, #,</>). (b) Surfaces dSi, dS 2 and dSs subtend the same solid angle dQ 
at point P. (c) Surfaces dS 2 and dS± are parts of the closed surface of a tetrahedron, with d5 2 
subtending solid angle d£l at point P, whereas dS± subtends solid angle —dfl at point P. 


Another way of considering the electric flux through a spherical surface is in terms of the solid 
angle subtended by surface element dS at the centre of the sphere. The solid angle can be 
defined in an analogous way to an ordinary angle. An arc of length £ which is part of a circle 
or radius r subtends angle 6 — i/r radians (rad) at the centre, so that if the arc were closed 
(i.e. a circle) the angle would be 6 — 2nr/r — 2n rad. In a similar way area dS which is part 
of a spherical surface of radius r subtends solid angle dfl = dS/r 2 steradians (sr) at the centre, 
so that if the surface were closed (i.e. a sphere) the solid angle would be Q = 47 rr 2 /r 2 = 47 t sr. 
Hence the solid angle subtended at the centre by surface element dS in Fig. 1.3(a) is 


dVt 


dS 


sin 6 dd d(j) 


|d(cos 0 )| d(j). 


(1.6) 


In Fig. 1.3(b) the solid angles subtended at point P by small surfaces dS\, dS 2 and dS% are 
obviously the same - it is the solid angle dVt within the cone which has its apex at P. Notice 
that only dS\ could be part of a spherical surface centred on P, in this case having radius i?i, 
whereas dS 2 and dS% are not centred on P. However, what counts is the projected area, so that 


dQ i 


dSi • Ri 
— 


gK4 2 


dS 2 • R 2 


3 


dSs • R3 

— r! - 


(1.7) 


where Ri, R 2 and R 3 are separation vectors from P to dS 1 , dS 2 and dS 3 . For dSi, d5 2 and 
dS 3 in Fig. 1.3(b) dVt\ = dfi 2 = dQ 3 = dfl. 
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Supposing dSi , (IS 2 and dS 3 are parts of closed surfaces Si, S 2 and S 3 of arbitrary shape 
surrounding P, then the sold angles these surfaces subtend at P are 


dSi • Ri 


Si 


R\ 


dSo • R2 


S 2 


R\ 


dS 3 • R3 


5 3 


Rl 


— 0 dfl — 4:7V 


(1.8) 


On the other hand, if point P isn’t inside the closed surface (e.g. P is outside the tetrahedron 

f dS • R 

in Fig. 1.3c) then 0 ——=— = 0; notice that the contributions from dS 2 and dS 4 cancel each 

Js R 

other out as dQ 4 = —dQ 2 . 


Returning to Fig. 1.3(b), the electric fluxes due to charge q at P out through surrounding closed 
surfaces Si, S 2 and S 3 (of which dS\, dS 2 and dS 3 are surface elements) are 

Q / dS 1 • R _ q f dS 2 • R _ g f dS 3 • R _ q / _ q , . 

4tt6:o /s', 47r ^o Js 2 Rl 4:7 T£q J Ss Rl 47T£ 0 / £0 ’ 

and this is the same for all closed surfaces surrounding charge q. 

We can also think of this in terms of the number of field lines leaving a charge per unit solid 
angle being isotropic (the same in all directions), and the electric flux through a suface being 
proportional to the number of field lines which in turn is proportional to the charge. Then the 
electric flux out through a closed surface is proportional to the solid angle it subtends, which 
is 47 r if the charge is inside the closed surface, or zero if the charge is outside. If we double the 
charge the flux would double, or if we add additional positive charges inside the closed surface 
the electric flux would always equal the sum of the charges inside, divided by £q. 

Since the direction along a field line is identical to that of the field direction, field lines must 
start on positive charges. Similarly, they must end on negative charges (or extent to infinity). 
Examples of field lines for different charge configurations were already shown in Figs. 1.2(d)-(f). 
Consider closed surface S surrounding all of the charge in Figs. 1.2(d)-(f). For Fig. 1.2(d) the 

net number of field lines leaving S is zero as is the net charge inside, so — (q — q)/so = 0. 

For Fig. 1.2(e) — (q + q)/s 0 = 2g/^o, while for Fig. 1.2(f) for which twice as many field lines 

leave surface S as enter it <$>e = (2 q — q)/£o = q/£o- 
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1.1.2 Gauss' law 


This law, named after German mathematician and physicist Johann Carl Friedrich Gauss 
(1777-1855), is a formalisation of the result of the previous section. Gauss’ law in integral 
form is 


(D E • dS — Qencl 

Js £o 


( 1 . 10 ) 


where Q enc i is the net charge enclosed by surface S. Applying Gauss’ theorem (Eq. C.19) we 
get 



V • E d 3 r — 



( 1 . 11 ) 


where surface S encloses volume V, and p( r) is the charge density (Cm 3 ). Equating integrands 
we obtain Gauss Law in differential form 

V • E = —. (1.12) 

£o 
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Where the symmetry permits, Gauss’ law in integral form (Eq. 1.10) gives a powerful method of 
finding the electric field. The trick is to find a volume satisfying: (i) surfaces are either normal 
to or parallel to the electric field, (ii) where the electric field is normal it is constant over that 
surface, (iii) the charge inside can be easily calculated. In the case of surface charges it is useful 
to construct a Gaussian pillbox with narrow sides and top and bottom surfaces parallel to the 
surface containing the charge. Classic examples include the infinite plane with uniform surface 
charge density a (C m -2 ) shown in Fig. 1.4(a), and the line charge density A (C m -1 ) of infinite 
length shown in Fig. 1.4(b). 


(a) infinite plane with 


(b) 




Figure 1.4: (a) Infinite plane with surface charge density cr, showing a Gaussian pillbox strad¬ 
dling the plane and whose upper and lower surfaces are parallel to the plane, (b) Section of 
infinite line charge surrounded by a Gaussian cylindrical surface. 

For the case of the infinite plane, from symmetry arguments the direction of the electric field 
must be normal to the plane, and for positive charge, E must be pointing away from the plane. 
Then, for area S of the infinite plane the charge inside the pill-box is aS so that 


|E|25 = —. 
£o 


E= ^ s ' 


(1.13) 


where the 2 S comes from the pillbox having both an upper and a lower surface. 

For an infinite line charge A we can use Gauss’ law for a cylindrical Gaussian surface of cylindrical 
coordinate radius p (Fig. 1.4b). Length L of the cylinder contains charge XL so that 


|E|27rpL —. 


E = 


A 


27 Tp£ 0 


P • 


(1.14) 


1.1.3 Conductors and insulators 

To understand why some materials are conductors and some are insulators, requires a brief 
excursion into quantum mechanics. In a solid, it is the outermost electrons in each atom, the 
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valence electrons, which are potentially free to move in the quantum well of the volume of the 
solid which has embedded in it a 3D periodic potential due to the ions in the crystal lattice. 
As a result, the discrete energy levels of the atoms in the solid merge into separate bands of 
closely spaced energy levels, possibly with a substantial energy gap between bands. 


The energy levels are populated with all the electrons in the solid starting with the lowest levels, 
but with a maximum of two electrons per level (one with spin up and the other spin down) as 
allowed by the Pauli exclusion principle. As one level is filled, remaining electrons go in the 
next level, and so on. The energy of the highest filled level is called the Fermi energy, Ep. Of 
course, at finite temperature there is a spread of electron energies above and below the Fermi 
energy, except when the band is completely filled and the energy band gap is too large. Where 
Ep is within a band determines whether the material is a conductor, semiconductor or insulator 
(Fig. 1.5). 


vacuum 


conduction 

band 


valence 

band 


✓ electron 
/ jumps 
/ band 
/ gap 



insulator semiconductor 



Figure 1.5: Electronic band structure of insulators, semiconductors and conductors. Red 
represents filled levels, and yellow represents vacant levels. 


Materials with only one valence electron per atom typically half fill the valence band, whereas 
materials with two valence electrons per atom fill the valence band. If the Fermi energy is in the 
middle of the band then lots of electrons with energies near the Fermi energy are able to move to 
nearby vacant levels in response to an electric field, as this takes only a tiny amount of energy, 
and are therefore able to move freely throughout the solid, making the material a conductor. 
Also, if there are exactly enough electrons to fill the valence band but the conduction band 
overlaps it, then electrons can easily move to the conduction band and the material is also a 
conductor. 


In contrast if the Fermi energy is at the top of a band, i.e. the band is full and there are no 
electrons in the next band, the material is an insulator. In insulators electrons are trapped in 
their energy levels; it takes a lot of energy (the band gap energy) to jump to the next level. 
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In semiconductors the band gap is sufficiently small that a few electrons are able to jump 
the band gap as a result of thermal excitation and/or in response to an electric field. In the 
process they leave behind vacant levels in the valence band or “holes” which act like positive 
electrons and contribute to conduction, just as the electrons that jumped to the conduction 
band, thereby causing the material to become semiconducting. Semiconductors used to make 
electronic components are typically “doped” with impurity ions which introduce energy levels 
within the band gap. 

In a conductor, electrons move in response to an electric field, which in electrostatics is due to 
other charges, and will keep moving until they reach the surface where they will redistribute 
themselves until the electric field inside the conductor vanishes. Thus, in electrostatics the 
entire conductor has the same potential, and its surface is an equipotential surface. 

Early experiments involved rubbing glass with silk, or rubbing amber with fur. What is going 
on here is that ruptured chemical bonds on the surface of one material leave a charge imbalance, 
and when that material is in contact with a different material with a different charge imbalance 
electrons will flow from one material to the other. It is the contact area that is important in 
charge transfer, rubbing just increases the contact area. The consequent displacement of 
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electrons from the glass to the silk leaves the glass positively charged, and from the fur to 
the amber leaves the amber negatively charged. Touching the charged glass to an isolated 
conductor would cause some electrons to flow from the conductor to the glass leaving the 
conductor positively charged, whereas touching the conductor with charged amber would cause 
some electrons to flow from the amber to the conductor leaving it negatively charged. 


1.2 The electrostatic potential 


The work done, W, in moving a particle of charge Q from point A to point B along an arbitrary 
path in the electric field of charge q located at r q is 


W = - [ QE(r) • dr 
Ja 



qH 

4t ts 0 R 2 


• dr 




Q 


4:7T£qR 


•dr (1.15) 


where R = (r — r g ), and we have used an identity for Vi? 1 (see Appendix B.l). Because 
(Vi?" 1 ) • dr is an exact differential, W is independent of the path taken from A to B. Ap¬ 
plying the principle of superposition for electric fields, we see that this would be true for any 
distribution of charges, and so for any electrostatic field. 


We have just shown that E is conservative in electrostatics. A conservative field is one for 
which its curl is zero, or equivalently it’s line-integral around a closed loop is zero. Thus, we 
may define an electrostatic potential difference 

w r B 

Vb-V a = — = - E(r) • dr 
V Ja 

In general, we may define the electrostatic potential 

V(r) = - f E(r / ) • dr' (1.17) 

J r 0 

such that 


(1.16) 


E(r) = -W(r). (1.18) 

The choice of the reference point ro, where the potential is defined to be zero, is arbitrary and 
does not effect the value of E(r). It is often chosen to be “at infinity”. If we set ro = oo we get 
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the potential for point charge q at the origin, and for point charge q' at location r' 

(i- 

respectively, where, as usual, R = r — r'. Applying the principle of superposition, the potentials 
due to volume and surface charge distributions p( r) and cr(r) are 


V(r) = 


47T£o r ’ 


and V (r) — 


q 


4:7T£()R 


V(r) 


1 

47T£o 


p( r 0 

R 



R 


dS'. 


( 1 . 20 ) 


Imagine plotting equipotential surfaces - if V only depends on say x and y you could plot a 
contour map and think of high values of V corresponding to hills, and low values as valleys. 
Then E at a particular location points down the steepest gradient, and is perpendicular to the 
contour lines. This also means that the electric field (and field lines) are normal to equipotential 
surfaces. 


1.2.1 Energy stored in the electric field 

Consider charge, q\ located at iq. The work done to move a second charge q 2 from infinity to 
r 2 , where the potential (due to the other charge) is V 2 = F(r 2 ), is 


W — q2^ / 2 — Q2 


Qi 


47r£ 0 |r 2 - ri| 




In general the work done to bring together a group of N charges is 


1 N 
w = 


( 1 . 22 ) 


2=1 


where Vi is the potential due to all charges except charge By analogy, the work done to form 
a continuous charge distribution is 


W 


= l[ p(r)V{r)d 3 r. 

^ J all space 


(1.23) 
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This formula cannot be used for point charges as the potential would be infinite at the locations 
of the charges. In nature this is not a problem as point charges (e.g. electrons) become fuzzy 
due the wave nature of particles (Heisenberg uncertainly principle). 


energy is 

stored in the electric field, as 

we can find by rewriting Eq. 1.23 as 


W = \ 

f £ 0 V • E(r) V(r) d 3 r, 


(1.24) 


J all space 



1 

o 

[ e 0 V • [E(r)V(r )}d 3 r - 

- \ [ eoE(r) • W(r) d 3 r, 

(1.25) 

Z 

J all space 

^ -'allspace 


1 

" o 

/ £ 0 [E(r)F(r)] • dS + \ 

[ e 0 E(r) • E(r) d 3 r, 

(1.26) 

z 

J S at oo ^ 

J all space 


- / 

$E(r) • E(r) d 3 r , 


(1.27) 

./allspace z 




where we have used a product rule (Eq. 1.8) to integrate by parts, and Gauss’ theorem (Eq. C.19) 
to replace the volume integral of V • (EE) over all space by the surface integral of EE over the 
closed surface at infinity — this is zero as EV ~ r~ 3 while S ~ r 2 as r oo. Hence, we see 
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that the energy density stored in an electric field is 


«E(r) = |E(r).E(r) = | (E(r)] 2 . 


(1.28) 


Summary of Chapter 1 

Coulomb’s law 

— Force on charge Q by charge q a distance R away is proportional to Qq/R 2 . 

— Force on Q is due to the electric field E of q acting on Q, i.e. F = QE(rg). 

— Electric field is defined as the force per unit test charge (unit: V m _1 ). 

— Electric field obeys the superposition principle. 

— For a point source q (C) or volume charge density p(r') (C m -3 ) at the source point r' the 
electric field at the field point r is 


E(r) = 


gR 


4t t£ 0 R 2 ’ 


E(r) ^/ 


P(r'> 


where R = (r — r 7 ). 

— £o = 8.85418782 x 10 -12 F m _1 is the permittivity of free space. 

Electric flux and Gauss’ law 


(1.29) 


— Electric flux through surface S is defined by — / E • dS. 

Js 

— Gauss’ law: (f E • dS — —Q e neb V • E = —. 

Js £ o £ o 


Electrostatic potential 


E is conservative in electrostatics and we may define an electrostatic potential 


V(t) = - f E(r 7 ) • dr 7 , E(r) = -VE(r). 

j 1*0 


(1.30) 


23 


Download free eBooks at bookboon.com 





Essential Electromagnetism 


Electrostatics 


— The choice of reference point 1 * 0 , where the potential is defined to be zero, is arbitrary and 
does not effect E. It is often chosen to be “at infinity” 


Energy stored in the electric field 


— The work done to bring together N charges, or a continuous charge distribution p(r) is 


1 N 1 r 

^ = w = o\ 

^ ^ J all space 


p( r)V (r) d 3 r. 


(1.31) 


The energy density stored in an electric field is ue(?) = ^E(r) • E(r) 


Exercises on Chapter 1 


1-1 The surface of a non-conducting sphere of radius a centred on the origin has surface charge 
density cr(a, #,</>) = crocos# and is uniformly filled with charge of density po- Find the 
electric field at the origin. 

1-2 A spherically symmetric charge distribution has the following charge density profile 


p{r,0,4>) 


Po (r < a ) 

/9 0 (r/a) _/3 (r > a) 


(1.32) 


where ft is a constant (2 < < 3). Find the electric field and electrostatic potential 

everywhere. 

1-3 The electric field is given by E(r) = Eq cos(z/zo) exp(— r / ro ) r , where zq and ro are con¬ 
stants. Find the charge density. 

1-4 If we had a point charge q at the origin we might choose the reference point to be some 
point at an arbitrary distance ro (usually infinity) from the origin. Then if we wish to 
find V(r,0,(j)) it would be convenient to have the reference point at ro = (ro, #,</>). Al¬ 
though obtaining the potential in this case is trivial, and one would usually just write it 
down, obtain the potential by carrying out explicitly the line integral for an appropriately 
parameterised curve. 

1-5 The electric field is given by E(r) = Eq cos(z/zo) exp(— r / ro ) r , where zq and ro are con¬ 
stants. Check whether or not the electric field is conservative. If it is conservative find 
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the potential, if it isn’t suggest how it may be possible to find the electrostatic part of the 
electric field (if present) and the corresponding electrostatic potential V(r). 


1-6 


How much work must be done to assemble: (a) a physical dipole made of charge +q and 
charge — q separated by distance d, (b) a physical quadrupole made up of four charges +q, 
—q, +q and —q on successive corners of a square of side d, and (c) a physical quadrupole 
made up of four charges —q, +q, +q and — q equally spaced apart by distance d on a 
straight line (see diagram below). 


(a) 


# + v ' ^ ( 

l -— ii 


+ 

d 



+ — 

ii - < 

>- o 


1 


3 


2 


(c) 


d. 


- + + 

#-— •- * 

1 2 3 


4 


1-7 (a) Use Gauss’ law in integral form to find the electric field due to charge density p(r) = 
po exp(—r/ro), and (b) check that you obtain the original charge density by taking the 
divergence of the electric field you find. 

1-8 An isolated conducting sphere of radius a has net charge Q. Find how much work was 
done to charge the sphere using two different methods: (a) from the charge on the sphere 
and its potential, (b) by finding the energy stored in the electric field. 
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2 Poisson's and Laplace's equations 

Learning objectives 

— To be able to derive Poisson’s equation from Gauss’ law, and know that Laplace’s equation 
is a special case of Poisson’s equation which applies in charge-free regions. 

— To know that the method of images for solving Poisson’s equation works by substitut¬ 
ing appropriate image charges for conducting boundaries, and to be able to obtain the 
potential for simple cases where there is plane, spherical and cylindrical symmetry. 

— To be able to solve Laplace’s equation in Cartesian coordinates in 3D for the potential 
inside a rectangular box where the potential or surface charge density is specified on all 
surfaces, and also for the case where the potential only depends on two of the coordinates. 

— To be able to solve Laplace’s equation in spherical coordinates inside and outside a spher¬ 
ical surface where the potential or surface charge density is specified, using the Legendre 
polynomials in the case of no azimuthal dependence of the boundary conditions, and using 
spherical harmonics for the general case. 

— To learn how to use a finite-difference numerical method for solving Laplace’s equation. 


2.1 Poisson's equation 


Named after Simon-Denis Poisson (1781-1840), Poisson’s equation is obtained from Gauss’ law 
by taking the divergence of E = -W, 

V-(-W) = (2.1) 

£0 

giving Poisson’s equation 

V 2 l/(r) = (2.2) 

£0 

We already know its solution for the boundary condition V(r —» oo) = 0, i.e. Eq. 1.20 


V(r) 


1 

47T£o 



pV) 

r — r 



(2.3) 
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As an aside, we can use Poisson’s equation and the potential due to a point charge to check 
an identity involving the 3D Dirac delta function 5 3 (r) which is described in Appendix C.7. A 
point charge q located at r q is represented by charge density p(r) = qd 3 ( r — r q ). But for this 
case the potential is given by Eq. 1.19, so 


V 2 



.-. V 2 



~s 3 ( r-r,), 
£o 

-47r5 3 (r - r q ), 


(2.4) 

(2.5) 


in agreement with Eq. C.38. 


Often the boundary conditions differ from V(r oo,0,0) = 0, and we will need to solve 
Poisson’s equation or Laplace’s equation taking account of boundary conditions determined by 
the problem. The usual types of boundary condition are Dirichlet boundary conditions for which 
V is specified at a boundary (e.g. a system of conductors held at various known potentials), 
and Neumann boundary conditions for which dV/dn | boundary is specified at a boundary. Note 
that dV/dn represents the rate of increase of V with distance in the direction normal to the 
surface, and that E • n = — dV/dn. From Gauss’ law, this is equivalent to specifying the surface 
charge density at the boundary surface. In Appendix E it is shown that the solution of Poisson’s 
equation is unique if V is specified anywhere on the boundaries, or is unique up to an additive 
constant potential if E • n is specified on the boundaries. 


2.2 Method of images 

2.2.1 Point charge above a grounded plane conductor 


The method of images is a powerful way of solving Poisson’s equation where the symmetry of 
the problem permits. The simplest case is that of a charge —g, say, some distance away from a 
grounded plane conductor as shown in Fig. 2.1(a). If we have two equal but opposite charges in 
empty space then the potential in the plane midway between the charges will be zero (Fig. 2.1b), 
giving the same electric field above that plane as above the grounded conductor with the single 
charge —q. It is as if at the mirror-image point below the conducting plane there is another 
charge, the image charge , q^ which in the case a plane conductor is the same magnitude but 
opposite sign to the real charge. 
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(a) 



(b) 



Figure 2.1: (a) Point charge — q above a grounded plane conductor showing field lines (blue 
curves with arrows), equipotential contours (red curves) and the location of the image charge, 
(b) Field lines and equipotentials of two real charges — q and +q in the same locations as the 
real charge and image charge in part (a). 


As an example, consider the case of a grounded conductor in the xy plane with charge — q 
located at (0, 0, +zo). The potential and electric field for z > 0 can be calculated as the sum of 
those due to the real charge — q and image charge -\-q at (0, 0, —zq), 
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V(x, y, z > 0) 
E {x,y,z > 0) 


+ 


+q 


^TTSo^x 2 + y 2 + (z - z 0 ) 2 4tt£o^/x 2 + y 2 + (z + z 0 ) 2 ' 


( 2 . 6 ) 


-q xx. + yy + (z - z 0 )z +q x5i + yy+(z + zp )z 

47T£o [x 2 + y 2 + (z — ^ 0 ) 2 ] 3/2 47T£:o [x 2 + y 2 + (z + zq ) 2 ] 3 / 2 


The resulting field lines and equipotential surfaces (equipotentials for short) are plotted in 
Fig. 2.1(a). Notice that at the surface of the conductor the electric field lines (and E) are 
normal to the surface. The reason for this is that electric field lines start or end on charges, 
and that if E was not normal to the surface at some location, the charges on the surface 
would experience sideways forces redistributing them across the surface until E became normal 
everywhere on the surface. 


Surface charge will be present on the conducting surface, as we can easily show using Gauss’ law. 
Imagine a Gaussian pillbox similar to that in Fig. 1.4(a) but with the important difference that 
the plane is the surface of a conductor and the volume below the plane (where the lower part 
of the pillbox is) is inside the conductor. Because in electrostatics E = 0 inside a conductor, 
there is only a contribution to j> E • dS from the upper surface of the pillbox, and the surface 
charge density on a conductor must therefore be in general 

dV 

o-(r) = e 0 E • n = e 0 (-W) • n = ~ £ o~q^, (2-8) 

where dV/On is the directional derivative normal to the surface and is defined as (W) • n. 

Returning to the example above, at the surface z = 0, and so the electric field from Eq. 2.7 
becomes 

(x, y ,0) 2 7T£o [x 2 + y 2 + ^o]3/2 ’ 

and hence the surface charge density on the conductor is 

r n\ V z ° q zo _ , , 

= < 2 ' 10 > 

where s = ( x 2 + y 2 ) 1 / 2 . Notice that the induced surface charge density in this case is positive 
everywhere, and this is because the charge near the conductor which induced the surface charge 
is negative. 


(2.9) 
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We can integrate the surface charge density over the conductor to find the total induced surface 
charge, 


Qs 


-r 

J o 


2ns <r(s, zo)ds = — 
2n 


2ns 


zo 


0 (s 2 + Z%) 3 / 2 


(IR = —q 


zo 




+q. 

( 2 . 11 ) 


Finally, the force between charge — q and the grounded plane must be the same as would exist 
between the real charge — q and a real charge of magnitude equal to qi (in this case +q) at the 
image charge’s location 


i <m = _i <f 

47T£o R 2 qqi 47T£ 0 (2^o) 2 


( 2 . 12 ) 


where R qqi is the distance between the image charge and the real charge, in this case R qqi = 2zq. 


2.2.2 Point charge outside grounded spherical conductor 


The method of images may be used to find the potential due to a point charge q outside a 
grounded conducting sphere of radius a (as in Fig. 2.2a) by replacing the sphere by a fictitious 
image charge, qi, inside the conductor. From symmetry arguments, the image charge qi must 
be on the line from sphere’s center to charge q , and the sum of the potentials due to the real 
charge and image charge must be zero on the surface of the sphere, e.g. at point P, 


V(r) = 


i fq , Qi 


47T£o 

1 


\d b)’ 


+ 


qi 


47re ° y y/a 2 + z 2 - 2 ~azg cos 9 \fa? + - 2azj cos 9 J 

1 ( q/a qi/zi 


+ 


4tT£o l yf\ + (z 9 /a) 2 — 2 (Zq/a) COS 9 Y 


(2.13) 

(2.14) 

(2.15) 


By inspection, the boundary condition V = 0 on the surface of the sphere is satisfied for any 
6 -value by q/a = —qi/zi and z q /a = a/zi, from which we obtain the location and magnitude of 
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the image charge 


Zi = a 2 /z q i qi = —(a/zq)q. 


(2.16) 


The resulting field lines and equipotentials are plotted in Fig. 2.2(b). 



Figure 2.2: Point charge outside grounded spherical conductor showing (a) position of image 
charge, and (b) field lines and equipotentials. 
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2.2.3 Line charge outside grounded cylindrical conductor 


We previously used Gauss’ law to find E due to an infinite line charge (Fig. 1.4b, and Eq. 1.14) 


e = p- 

2npeo 


(2.17) 


From the electric field we find the potential 


V{r) = - f E(r') • dr' 
j r 0 


(2.18) 


where ro is the reference point where the potential is defined to be 0. Because of the cylindrical 
symmetry about the line charge, we can define the reference “point” to be anywhere on a 
concentric cylinder of arbitrary cylindrical radius po 5 and so the potential is 


V{r) = 


f 

J po 


X 


27 xpeo 


P ' 0 dp) = 


X f p dp 

2tt£o J P0 P 


2tt£ 0 


In I — 


(2.19) 



Figure 2.3: Line charge outside grounded cylindrical conductor showing position of image line 
charge, and geometry for finding the potential. 


Next consider the potential due to two parallel line charges Ai and A 2 . Suppose that for line 
charge 1 the reference “point” is arbitrarily set at cylindrical radius away from it, and 
for line charge 2 the reference “point” is arbitrarily set at cylindrical radius p q y away from 
it. Clearly the two reference “points” can not coincide, and so one should ask whether or not 
this matters, but as we shall see it doesn’t. Let’s just go ahead and apply the principle of 
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superposition and write down the potential due to the two line charges 


Ai 

m( p ° ) )+ A 2 in( 


( 2 . 20 ) 

2tt£o 

\ Pi J 2tt£q \ 


1 

2tt£o 

(Ai In pp + A 2 ln/4 2 ^ 

) 27r£o (Ailnpi+A 2 lnp2) 

( 2 . 21 ) 

Eq - 

(Ai In pi + A 2 lnp 2 ) • 

2ire 0 

( 2 . 22 ) 


where p\ is the distance of line charge 1 to the field point r, and p 2 is the distance of line 
charge 2 to the field point. Note that in Eq. 2.21 the first term involves only constants, the 
constant line charge densities Ai and A 2 and the arbitrary constants p^ and p^ \ and so may 
be written as an arbitrary constant potential Vo whose value will not affect the electric field or 
the potential difference between two points. 

Equipotential contours and field lines are plotted in Fig. 2.4(a) for the the case of two parallel 
line charges with equal but opposite line charge densities, +A and —A, for which 

= < 2 - 23 > 

Notice that the equipotential surfaces for this case are of circular cross section. This tells us 
immediately that for the case of a line charge and a parallel grounded conducting cylinder we 
can use the method of images with an equal but opposite image line charge inside the conducting 
cylinder and parallel to the real line charge. 

In Fig. 2.3 The line charge density A (C m -1 ) is parallel to the conducting cylinder’s axis and 
distance d from it. We have just seen that we can use an image line charge —A which from 
symmetry must be parallel to to and in the plane containing the cylinder axis and the real line 
charge A. One can show (Exercise 2—3) that the cylinder’s surface is an equipotential if 

di = a 2 /d. (2.24) 

Equipotential contours and electric field lines are shown in Fig. 2.4(b). The result above may 
also be applied to the case of two parallel conducting cylinders, e.g. a twin wire transmission 
line having equal but opposite charges (Fig. 2.4c) in order to work out the capacitance per unit 
length (Exercise 2-4); the capacitance of two conductors with potential difference V and having 
charge +q on one and — q on the other is C = q/V and will be discussed further in Chapter 4. 
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Figure 2.4: Electric field lines and equipotentials around: (a) two parallel line charges +A 
and —A, (b) one line charge +A and a parallel conducting cylinder, (c) two parallel conducting 
cylinders with equal and opposite charge. In (b) and (c) the dots show the positions of image 
charges. 
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2.3 Laplace's Equation 


Named after Pierre-Simon, Marquis de Laplace (1749-1827), Laplace’s equation is simply Pois¬ 
son’s equation for p(r) = 0, and in the three main coordinate systems is given by 


Cartesian: 


Cylindrical: 


d 2 V d 2 V d 2 V 
dx 2 + dy 2 + dz 2 

1 d ( dV\ 1 d 2 V d 2 V _ 
p dp \ dp ) + p 2 d(f) 2 + dz 2 


on . , 1 d ( 2 dV\ 1 

Spherical: V^g-lr— +^ n/)m 


d ( . n dV\ 
sin 9-—- I + 


1 d 2 V 


dO J r 2 sin 2 0 deft 2 


= 0. 


(2.25) 

(2.26) 

(2.27) 


It must be solved together with boundary conditions appropriate to the problem. 


2.3.1 3D Laplace's equation in Cartesian coordinates 

Using the method of separation of variables we will find the general solution for the potential 
inside a rectangular box where the potential is specified to be zero on all but one of its surfaces. 
The particular solution can then be obtained by applying the boundary condition that the 
potential must be that specified on the surface which has V ^ 0. Having solved for this case, 
we can do the same for V being non-zero on each of the other five sides, and then use the principle 
of superposition to get the full solution. The geometry to be used for solving Laplace’s equation 
is shown in Fig. 2.5(a) together with a cartoon (Fig. 2.5b) indicating how the full solution may 
be built up from six solutions, each being for when the potential is specified on one of the 
surfaces and is zero on the other five. 


(a) kZ (b) 



Figure 2.5: (a) Geometry for solving Laplace’s equation in 3D in Cartesian coordinates for the 
field inside a rectangular box where the potential is specified on the six surfaces, (b) A cartoon 
indicating how the full solution may be built up from six solutions. 

Using the method of separation of variables we try to find a solution of the form V(x,y,z) = 
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X(x)Y(y)Z(z ) and substitute this into Laplace’s equation, which becomes 


d 2 X(x) 
dx 2 


i + W + -&) xix)ny)z(z> 

Y(y)Z(z) + X(x)^Xz(z) + X(x)Y(y) d2Z{z) 


dy 2 


dz 2 

X"{x) Y^y) Z^z) 
X(x) Y(y) Z{z) 


= 0 , 

(2.28) 

= 0 , 

(2.29) 

= 0 . 

(2.30) 


Note that the 1st term depends only on x, the 2nd only on y and the 3rd only on z. The 
only way that the equation can hold as x, y or z are varied independently is if each term is a 
constant, and the three separation constants add to zero, 


X"(x) Y^y) Z"(z) 

X(x) '• Y(y) 2 ’ Z(z) 


S 3 , where Si + S 2 + S 3 = 0 . 


(2.31) 


At least one of the three separation constants must be positive and at least one must be negative. 


For our hollow rectangular box (Fig. 2.5) let’s choose the surface with non-zero potential to be 
that at z — c, on which the potential is specified to be V{x,y,c) — Vo(x,y). It makes sense to 
have Si and S 2 negative so that the equations in x and in y will have periodic solutions, making 
it easy to satisfy their boundary conditions. Then the three equations to solve are 


X"(x) = -a 2 X(x), Y"(x) = -/3 2 Y(y), Z"(x) = 7 2 Z(z), (2.32) 

where y 2 = a 2 + /3 2 . These equations have well-known solutions 

X(x) = a\ sin (ax) + a 2 cos(crx), (2.33) 

Y ( y) = 61 sin(/3y) + b 2 cos(/3y), (2.34) 

Z(z) = ci exp (72:) + c 2 exp (—72:), (2.35) 


as can easily be checked by substitution. For Eqs. 2.33 and 2.34 the boundary conditions, i.e. 
V (0, y,z) — V (a, y, z) = 0 and V (x, 0 ,z) — V (x, b, z) — 0, require a 2 = b 2 = 0 and that 

a = a k = ^~, ft = fa = y, and so j 2 = 7 kl = a 2 k + p 2 , (2.36) 

where k and l are integers. For Eq. 2.35 the boundary condition V(x, y, 0) = 0 requires c 2 = —ci, 
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and so Z(z) oc ( e iz — e lz )/2 — sinh( 7 z)\ actually, we could equally well have written the 
general solution of Eq. 2.35 in terms of the hyperbolic sine and cosine functions (Fig. 2.6), i.e. 


Z(z) = d\ sinh( 7 z) + g ?2 cosh( 7 z) = d\- 


alZ _ p~lZ 


+ G?2 


e iz + e ~iz 


(2.37) 


A useful property of the hyperbolic sine and cosine functions is that 


sinh(z) = cosh( 2 :) and — cosh (z) = sinh(^). 

Li/C Li/C 


(2.38) 


Then, the general solution to Laplace’s Equation for the case where V(x,y,c) 7 ^ 0 with the 
potential on all other surfaces being zero is 


v ( X,y,z ) = EE A k i sinh ( 7 k i z) sin (a k x ) sin ($ y). 


(2.39) 


k =1 1=1 



> Apply now 


REDEFINE YOUR FUTURE 

AXA GLOBAL GRADUATE 
PROGRAM 2015 


redefining /standards 
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x 


Figure 2.6: Hyperbolic sine and cosine functions compared with the positive and negative 
exponentials. 

The coefficients Am can be obtained using the boundary condition at z = c, i.e. by solving 

oo oo 

EE Aki sinh (p( k i c) sin (a k x) sin (/?; y) = V Q (x,y). (2.40) 

k =1 1=1 

The easiest way to solve this is by recognising that the left hand side of Eq. 2.40 is actually a 
2D Fourier sine series, i.e. 


V 0 (x,y) = sin 


k =l 1=1 


\ a J 


(2.41) 


for which the Fourier series coefficients are 


Bki — 


4 f a f b fk7Tx\ . (hry'. 

ab i l V ° (X ' y) “"I”) Sm (“ 1 dxdy - 


(2.42) 


Hence 


Aki = ~ 


B, 


ki 


n b 

Vo (. x , y) sin (a k x) sin ($ y) dx dy. (2.43) 

. 


sinh (7 m c) ab sinh (y ki c ) J 0 J 0 
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2.3.2 2D Laplace's Equation in Cartesian Coordinates 

If there is no dependence on z, we need only to solve Laplace’s equation in 2D for the region 
inside a rectangular cross-section pipe. Again the method of separation of variables can be used 
together with the principle of superposition to break down a problem with V ^ 0 on all four 
sides to 4 simpler problems with V ^ 0 only on one side. We need to decide after examining 
the problem whether it is sensible to have the periodic solution for X(x ) or for Y(y). Keep the 
negative sign for the equation for X ( x ) if the boundary conditions suggest it would be simpler 
to solve if we have a periodic solution in x , e.g. if V=0 at both boundaries in x. 

The general solutions are then of the form 

V(x, y ) = (A cos (kx) + B sin (kx)) [Ce ky + De~ ky ^j , (2.44) 

or equivalently 

V (x, y) — (A cos(kx) + B sin (kx)) [C' cosh (ky) + D f sinh(fcy)] . (2.45) 

Finally, the boundary conditions are applied, as for the 3D case, to find allowed values of k and 
the coefficients in the Fourier series for the potential on the side where V ^ 0, and so the full 
solution. For example, inside the rectangular region corresponding to 0 < x < a and 0 < y < b 
with K(x, b ) = Vo(x), and V—0 on the other 3 sides, 

oo 

V ( x , y) = A n sin (^L) sinh (^y) , (2.46) 

71=0 

and the boundary condition at y — b determines the coefficients A n in this Fourier sine series 

^ = l Va(x) sin (t 1 ) dx - <2 - 47> 


As an example we shall consider the case where a = b = 1 m, V(x, b) = 1 V and V = 0 on the 
other three sides. Then, 


n sinh (rar) 


/ sin ( mrx ) dx. 

J o 


n rm sinh {mt) 


(for odd n) 


(2.48) 
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Using Eqs. 2.46 and 2.48, and summing terms up to n = 19 we obtain the solution contoured 
in Fig. 2.7. 


1.0 


0.8 


0.6 

y 

0.4 


0.2 


0.0 

0.0 0.2 0.4 0.6 0.8 1.0 

x 

Figure 2.7: Contours of V for the the result using Eqs. 2.46 and 2.48, and summing terms up 
to n = 19. 
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We can obtain the electric field from the potential in the usual way 


Ex — 

dV _ 

- 4 E 

n, odd 

1 

■ cos(n7rx) sinh(n7r7/), 

dx 

sinh (mr) 

Ey = 

dV 

- 4 E 

n, odd 

1 

■ sin(n7rx) cosh(n7T7/). 

dy 

sinh (mr) 


(2.49) 

(2.50) 


Another example of solving Laplace’s equation in 2D is given in Exercise 2-5. 


2.4 Laplace's equation in spherical coordinates 

In spherical coordinates Laplace’s equation is 

J_d_ / 2 dE\ 1 d / . dV\ 1 d 2 V 

r 2 dr \ dr ) r 2 sin 9 dO \ dQ ) r 2 sin 2 6 dcj) 2 


(2.51) 


2.4.1 Solution for no azimuthal dependence 

Here we will first consider the case of azimuthal symmetry {dV/dcj) = 0) and look for solutions 
of the form V{r,Q) = R(r)Q(0). Then 


1 d 
R(r ) dr 


>dR(r) 

dr 


1 


d 


Q(9) sin 6 d6 


sin 9 


d©(0) 

dd 


(2.52) 


The left hand side depends only on r, and the right hand side depends only on 6 . Each side 
must equal the same separation constant C, and so we have two ordinary differential equations 


1 d 
R(r) dr 

1 d 


} dR(r) 

dr 


0(0) sin 6 dO 


sm9 ^r 


= C, 
= C. 


(2.53) 

(2.54) 


We will first tackle the radial equation (Eq. 2.53), and to simplify later working (we’ll see 
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shortly) set C = a(a + 1 ) and write it as 

d_ ^ = a(a + ( 2 . 55 ) 

This equation has general solution 

R(r ) = ,4r“ + £r _( “ +1) (2.56) 


as we can easily check by substitution. 


The angular equation (Eq. 2.54) can be re-written in terms of cos0, 


d 


d{ cos 6 ) 


(1 — cos 2 6 ) 


d& 


d(cos9) 


T cy(cy T 1)0 — 0, 


(2.57) 


which is known as Legendre’s equation. This equation can be solved by a power series method, 
and the series converges for — 1 < cos 9 < 1 for integer values of a = £ > 0. Since Legendre’s 
equation is 2 nd order, there are two linearly-independent solutions for each value of <a, of which 
one is infinite for cos 6 = ±1. The finite solutions are are called “Legendre polynomials”, 
P^(cos0), and the first 6 Legendre polynomials are 


Po(cos$) = 1, 

Pi (cos 6 ) = cos0, 

P 2 (cos$) = ^(3 cos 2 6 — 1), 

P 3 (cos 9) = - (5 cos 3 6 — 3 cos 9 ), 
z 

7^ (cos 0) = i (34 cos 4 9 — 30 cos 2 9 + 3), 

8 

P^(cos9) = -(63 cos 5 9 — 70 cos 3 9 + 15 cos 0), 
8 


(2.58) 

(2.59) 

(2.60) 
(2.61) 
(2.62) 
(2.63) 


and are plotted in Fig. 2.8. Notice that Pi(cos9) are normalised to 1 at cos 9 = 1, and that 
Pi(cos9) is an even function if £ is even, and an odd function if £ is odd. 


The Legendre polynomials form a complete orthogonal set 


J P(> (cos 9)P((cos 9)d cos 9 = j- 


Hl'i 


(2.64) 
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Figure 2.8: The first 6 Legendre polynomials. 


where Sij is the Kronecker delta which takes the value 1 if i 


j and 0 if i ^ j. Note the 
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similarity to the sine and cosine functions in Fourier series; this enables us to represent any 
function f(x) over the range -1 < r < +1 as a sum of Legendre polynomials 

op 4- 1 Z* 1 

f( x ) = y2Ce p e( x ), where C e = —-— / f(x)P^{x)dx. (2.65) 

r=o 2 J ~ l 

If f(x) is an even (odd) function, only terms with even (odd) £ will be non-zero. 

The General solution of Laplace’s equation with axial symmetry is then 

oo 

V(r, 9) = rt + B t r~ {e+1) )Pe( cos 9). ( 2 . 66 ) 

t =0 


The coefficients in the series are obtained from the boundary conditions. If the region inside 
the boundaries contains r = 0 then all B£ will be zero, but if the region inside the boundaries 
contains r — oo then (usually) all An will be zero. At the boundary (often a spherical surface 
with radius r — a) where the potential is specified, we can use the orthogonality of i^(cos$) 
to find the coefficients. Hence, for the region inside a spherical surface of radius r — a, where 
V(a,cos0) = Vo(cos$), we have 

op -u 1 f 1 

A e a e = Vo(cos 9)Pe(cos 9)d cos 9. (2.67) 

2 J -1 

Often it is simpler to re-write Vo (cos 6 ) as a series of Legendre polynomials and obtain Ai simply 
by comparing it with H(a,cos0) in Eq. 2.66. An example is given in Exercise 1-7. 


2.4.2 Solution including azimuthal dependence 


In this case the potential depends on both 6 and 0 at a spherical boundary. We will represent 
the angular dependence in terms of spherical harmonics 


YimM) 


'(2£ + l) (l- 1)! 

47 t {f + 1 )! 


pp (cos 9)e im(t> , 


( 2 . 68 ) 
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where 


d m 


m*) = - x 2 r/ 2 —p e ( X ) 


dx r 


(2.69) 


are associated Legendre functions. The first 9 spherical harmonics are given in Eq. 2.70, and 
the first 16 spherical harmonics are contoured in Fig. 2.9. 



Yi, o = \ —COS e, Yi,±i = =FW — sin6e ±l<t> . 


47T 


8n 


15 


Y 2fl = \/ —(3 cos 2 6-1), Y 2) ±i = T\I^cos0 sinfle^, y 2) ±2 


8 tt 


(2.70) 


]/ 


15 
327r 


sin 


Q e ±2i<P_ 



Figure 2.9: The first 16 spherical harmonics. Note: these polar plots are not equal-area 
projections. 

The spherical harmonics form a complete orthonormal set 


/*27T P 7T 

/ #/ dd sm6Y e * m (e,cf>)Y e ^(e,(f)) = 5e t ef5 mm f, (2.71) 

J 0 J 0 

and since we are integrating over all 6 and (j) we are actually integrating over all solid angles, 
i.e. 


dnY e * m (dA)Ye, r 


'{6 ,< t >) = 


(2.72) 
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Using spherical harmonics, the general solution of Laplace’s equation in spherical coordinates 
where there is azimuthal dependence is a straightforward extension of the case with no depen¬ 
dence on </>, 


(X) £ 

V(r,6,4>) = EE (At, m r e + B e , m </>). (2.73) 

£=0 m=—£ 

Suppose, then that there is a spherical surface of radius a centred on the origin on which the 
potential is specified as V (a, 0, </>) = Vo(0, </>). If we wish to find the potential outside the sphere, 
to keep V finite all A^ m must be zero in Eq. 2.73, and using Eq. 2.72 together with the boundary 
condition we obtain 

oo £ 

V(r > a,0,<j>) = E E r-V +1) Y ttm (0,<l>), (2.74) 

1=0 mn=—£ 
oo £' 

V 0 (8,(t>) = J2 E (2.75) 

£'=0 m'=—£' 
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h^ m (0,</>)F o (M) = E E (2-76) 

£'=0 m'=-£ f J 


/ oo f 

<mY e * m (o,<i>)Vo{0,<j>) = '52 e 

£'=0 m '=—£' 


Be,m = a (m) <t dn Y e * m (8,4>)V 0 (e, </>). 


Similarly, inside the sphere 


oo £ 


V(r < a, M) = E E ^)> 

^=0 m =—£ 

A e>m = a~ e <f dQYfte, cf>)V 0 (8, cf>). 


(2.77) 

(2.78) 


(2.79) 

(2.80) 


As an example suppose the potential on the surface of a sphere is 


V(a, 6,4 >) = Vi sin 0(1 + 2cos0) sin0, 


(2.81) 


and we need to find the potential outside. This relatively simple boundary condition can be 
re-written in terms of spherical harmonics as follows 


H(a,0,0) = Vi 


\[^ PMM) 


^(0, </))] +2 A /— [r 2 ,i(M) 


[>i>i(M)-vVi(0^)] 


(2.82) 


Hence, 


B 1,1 = Hi 




and £> 2 ,i — 


-Hi 
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giving 


a 2 /27r 


r* 


a 6 /27r 


v(r,e,4>) = — \Y ltl ( 6 ,d>) - Fi,-1 (e ,</>)] + 2v^J — [Fi 21 ( e,</>) - y 2( _i (t9 ,</>)] , 


a 2 a 3 

= Vi sin 6 cos (j) + 2Vi sin 0 cos 0 cos (f>. 


(2.83) 


2.5 Finite-difference method for Laplace's equation 


Numerical finite-difference methods may be used when analytic methods are impractical. They 
can be motivated by the equation of heat conduction in solids, which becomes Laplace’s equation 
for heat-source free regions in equilibrium. 

The equation of heat conduction in solids in a region where there are no heat sources is 
1 r)T 

V 2 r = —— (2.84) 

a at 

where T(r, t) is the temperature and a is the thermal diffusivity. If the boundaries are held at 
fixed temperatures then for any initial temperature distribution the temperature distribution 
evolves to a steady state at t —>> oo at which time 

V 2 T = 0, (2.85) 

which is Laplace’s equation. 

By analogy, in electrostatics we can solve Laplace’s equation by numerically “evolving” an initial 
potential in “time” to t —» oo using 

^ = V 2 v, (2.86) 

or, in practice, until V (r) stops changing. The potential will then have “relaxed” to a steady- 
state and obey Laplace’s equation. This equation does not represent the time-dependence of 
V (r) — we are still discussing electro statics ! 
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We shall solve Eq. 2.86 numerically by taking “time” steps St 
V(r ,t + 5t) = V (r, t) + [vV(r ,t)]5t, 


(2.87) 


calculating V (r, t) only at points on a 3D rectangular grid with spacing L and a “time” grid with 
spacing St. We shall obtain V 2 E by numerical differentiation, using the shorthand notation 

V^ k = V{x = iL,y = jL , z — kL, t — nSt ), (2.88) 

and carrying out numerical differentiation with respect to x 


dV 

dx 


(x+L/2,y,z,t) 

d 2 V, 


yw - is* 

L ' dx 


(x—L/2,y,z,t) 


- v^k 

L 


dx 2 


(x,y,z,t) 


V i+i,j,k + ,j, k ~ 2V i?j,h 

L 2 12 


Hence, we obtain the Laplacian of V in 2D and in 3D, 


(2.89) 

(2.90) 


7 D w 
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(2.91) 

6V&*) • 


Then a suitable iteration scheme in 2D would be 



(2.92) 

(2.93) 


The “time step” St is arbitrary, and so we can choose whatever we like for it. The method is 


simplified by setting St = L 2 /4 (for 2D) or St = L 2 /6 (for 3D), and so in 2D and 3D problems 
the iteration procedure is simply 



(2.94) 


(2.95) 


Thus, for each iteration we replace the potential on each non-boundary grid point by the average 
of the potentials at its nearest neighbours. 

As an example, we solve the 2D problem of the potential inside a pipe of square cross section 
(1 m by 1 m) where the potential is specified on each of the four sides: V (0, y, z) = 0, V (1, y, z)— 0, 
V(x,0,z)=0 and V(x, l,z)=l V. This is the same problem we used as an example for analytic 
solution (Fig. 2.7). We choose a grid spacing L=0.01 m, i.e. 100x100 cells (101x101 grid 
points), and after each iteration reset the outer grid points to the potentials on the boundary 
surfaces. At the two corners where V changes discontinuously from V=0 to V—l V we set the 
potential to be the average value, i.e. 1/(0,1, z)=V(l, 1, z)=0.5 V. The result is contoured in 
Fig. 2.10 after 100, 1000, and 10,000 iterations. 

Since the information on the potential can propagate from one boundary to another at a rate no 
faster than one cell per iteration, the number of iterations required to approach a steady-state 
must be much larger than the maximum number of grid points in any coordinate direction (101 
in this example). In the present example we see that the number of iterations required, ^10,000, 
exceeds this by a factor of 100. Because of the large number of iterations required, very efficient 
codes are needed, and should use array manipulation methods available in modern programming 
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100x100 cells, 100 iterations 100x100 cells, 1000 iterations 100x100 cells, 10000 iterations 



0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

XXX 

Figure 2.10: Numerical solution of Laplace’s equation in 2D for the boundary condition being 

V = 0 on all sides of the square except for the side at y = 1 where V = 1. Contour levels are 

V = 0.1,0.2,..., 0.9. 

languages, rather than using do loops, when computing the average potential around each grid 
point, and when resetting the boundary potential. The code used in this example was written 
in IDL and is given in Appendix G. 


Summary of Chapter 2 

Poisson’s and Laplace’s equations 

— Can be solved numerically, by the method of images in some cases, or analytically using 
Fourier methods if boundary conditions permit. 


Laplace’s equation in Cartesian coordinates 


/ d 2 d 2 d 2 \ 

\dx 2 dy 2 dz 2 ) 


V(x,y,z) 


= 0 . 


(2.96) 


— In 2D, inside rectangular region corresponding to 0 < x < a and 0 < y < b with V (x, b ) = 
Vo(x), and V—0 on the other 3 sides, 


o° 

V (. x, y) = A n sin (^cc) sinh (^y) 


n =0 


(2.97) 
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— The boundary condition at y — b determines the coefficients A n in this Fourier sine series 

2 1 /*® / r fij[ \ 

A " = «Sih(f7 l Vi(x)sin ( t x ) ix - (2 ' 98) 

— In 3D, inside rectangular region corresponding to 0 < £ < a, 0 <y <b and 0 < z <c with 
V(x,y,c) = Vo(x,y), and V—0 on the other 5 sides 


v ( X,y,z ) = EE Aki sinh (7 k i z) sin (a k x) sin ($ y ), (2.99) 

k =l 1=1 

rCHT 1 7T 

= —, A = y and 7^ = o| + /3f. (2.100) 

— The boundary condition at z = c determines the coefficients A k i in this 2 D Fourier sine 
series 


Aki 


4 

a b sinh ( 7 ^ c) 



sin (07 x) sin (/?/ y ) dxdy. 


( 2 . 101 ) 
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Spherical coordinates with no dependence on <fi 

1 d ( 2 dV\ 1 d ( . dV\ 

7>~dr V !fr) + r 2 sin 6 06 \°1)6 ) ~ ' 


( 2 . 102 ) 


— The general solution involves powers of r and Legendre polynomials P^(cos$) 


V(r,9) = J^(A*r* + B e r-( i+ V)P e (cosO). 
£=0 


(2.103) 


— If the space considered contains r — 0 then all Bn = 0, or if the space considered contains 
r — oo then all A£ = 0. 

— The boundary conditions at r — a determine the coefficients A,£ or B£ using orthogonality 
of Pi (cos $) 



cos 6) Pi (cos 9) d cos 9 


1 




(2.104) 


Spherical coordinates with dependence on <fi 


— The solution involves powers of r and spherical harmonics Yg m (0, <j>) 


oo £ 

V(r,0,<j>) = Y Y + 

£=0 m ——£ 


(2.105) 


— If the space considered contains r = 0 then all B^ m = 0, or if the space considered contains 
r — oo then all A^ m = 0. 

— The boundary conditions at r — a determine the coefficients A^ m or B^ m using orthogo- 
nality of Y e>m (6, 0) 


dQY e * m (e,4 ) )Ye, m/ (d 1 cf>) = 5e,ed n 


(2.106) 
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Finite-difference method 


The potential is defined on a rectangular grid with potentials on the boundary grid points 
fixed by the boundary conditions. 

— The potentials at the non-boundary grid points are set to reasonable values, e.g. average 
of boundary potential. 

— The potentials on the non-boundary grid points are replaced by the average of the poten¬ 
tials over the nearest 6 grid points (3D) or the nearest 4 grid points (2D). 

— The procedure immediately above is repeated until the potential stops changing. 


Exercises on Chapter 2 

2-1 Charge +q is located on the z axis a distance d/2 from a grounded plane conductor in 
the x-y plane. Find how much work was done to bring the charge to its current location 
using two different approaches: (a) the work done against the electrostatic force if the 
image charge were real and there was no grounded conductor, (b) the work done against 
the electrostatic force due to the induced surface charge 

/ _ -q z 

Vi d) q / 2 i 2 i 2\3/2 

2n (; x z + y Z + Z z ) 6 l z 

where z is the height of the charge above the plane. 

2-2 Charge +q is brought near to two orthogonal grounded conducting planes, one correspond¬ 
ing to the x-z plane and the other to the y-z plane. The charge is located at (a, 6,0). 
Find the work done in bringing the charge from infinity to its current location (a) by using 
the method of images to find the potential at the location of the real charge, and (b) by 
considering the force on the charge as it is brought from infinity. 

2-3 Show that the potential outside a long conducting cylinder of radius a in the presence of 
a long parallel line charge +A at distance d is identical to the potential of the line charge 
and a parallel image line charge —A at distance d{ from the cylinder’s axis towards the real 
line charge (see diagram below). [Hint: draw lines to point P from the two line charges. 
Use the cosine rule of triangles to write the two distances in terms of a, d and </> and use 
the formula the for potential due to a line charge, and superposition, to write a formula 
for the potential at P. Finally require that V does not change if (j) changes.] 


(2.107) 
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2-4 Find the capacitance of a two-wire transmission line comprising two identical parallel cylin¬ 
drical conductors of radius a whose axes are separated by distance D (see diagram below). 
You may use the result for the potential due to a line charge near a single cylindrical 
conductor to find the potential difference by replacing the cylinders by equal but opposite 
image line charges, +A and —A (C m _1 ). The capacitance of two conductors with potential 
difference V and having charge +q on one and —q on the other is C — q/V 



2-5 A region of space is bounded by three plane conductors as illustrated. Find the potential 
everywhere between the conductors. 


V(0,y)= 



2-6 Find the potential inside the rectangular region, 0 < £ < a, 0 < y < b and 0 < z <c with 
V(x,y,c) — Vo(x,y), and Y=0 on the other 5 sides, where 


Vq( x,y) = Vi sin 



sin 



(2.108) 
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2-7 The potential on a non-conducting sphere of radius a is given by 
V = Vo(3cos 2 9 + cos 9 — 1). 

(a) Find the potential and electric field inside the sphere. 

(b) Find the potential and electric field outside the sphere. 

(c) Find the surface charge density on the sphere as a function of 9. 

2-8 Consider a point charge on the z -axis at z = r'. Find l/(r, #,</>) in terms of Legendre 
polynomials for r > r '. 

2-9 The potential on the surface of a sphere is 

V{a,0,(j)) — V\ sin 9 sin </> + V 2 sin 9 cos 9 sin </>. (2.109) 

Find the potential inside the sphere. 
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3 Multipole expansion for localised charge distribution 

Learning objectives 

— To understand how the potential for a localised charge distribution may be approximated 
at large distances by the first one or two non-zero terms of the multipole expansion of the 
potential. 

— To be able to obtain the the monopole moment, dipole moment and quadrupole moment 
tensor in Cartesian coordinates for simple cases, and to be able approximate the potential 
using these. 

— To be able to find the electric field and potential due to electric dipoles. 

— To be able to perform a multipole expansion of the potential in spherical coordinates. 

3.1 Point charge on 2 : axis 


Consider the potential due to a point charge q on the z axis a distance r' from the origin 
(Fig. 3.1a). We can obtain the expression for the potential E(r) for r > r' as a solution of 
Laplace’s equation in spherical coordinates with azimuthal symmetry. For this problem we can 
use Coulomb’s law for the potential due to point charge as the boundary condition (Exercise 2- 
8). The result is found to be 

oo 

V(r,0,<(>) = ■— y^(r / )V~ ( ^ +1) P < g( cos6>). (3.1) 

71-60 £=0 


In Eq. 3.1, 9 is just the angle between r and r', so that cos 9 = r • ?' and 

oo 

= |r-E( r, ) V_( ' +1) ^( f - ?/ )- (3.2) 

71 "^° £=0 

Written this way the result is applies equally whether or not the charge is located on the z -axis 
(Fig. 3.1b). 

The result given in Eq. 3.2 may be extended immediately to a localised charge distribution 
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Figure 3.1: Geometry for potential due to: (a) point charge q 
away from the z axis, (c) arbitrary charge density p( r). 



(Fig. 3.1c) for which we will obtain the multipole expansion of the potential. For charge density 
p( r), the contribution of volume d?r' to the potential at r > r' is 


dV{ r) 


p(r')cPr' 

4:7T£q 


• r'). 

i=o 


(3.3) 


For a localised charge density p( r) the potential outside charge region is therefore 
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i 00 r 

y(r) = - y / p(r')(r'YPi(r -r')d 3 r'. 

47r£ ° U J 


(3.4) 


Far away from the charge region we can approximate this by first few terms 

y(r) « -- [ p(r')d 3 r' + --~ [ p(r / ) r (r • r') d 3 r' 

v y 4 t T£ 0 r J 4t T£ 0 r 2 J 

+ 4^3 /Kr')(r') 2 5 [3(f-f') 2 - 1] <iV + -.. (3.5) 

The first three terms are the monopole, dipole and quadrupole terms of the multipole expansion , 
and they are proportional to r _1 , r~ 2 and r -3 , with each having its own multipole moment. 
The lowest non-zero multipole moment is independent of the choice of coordinate origin, unlike 
higher multipole moments. 


3.2 The monopole moment 


The monopole potential is 


^mon(r) - 47r£or Jp( r ')d 3 r' 


Attest 


where the monopole moment is just the net charge q. 


(3.6) 


3.3 The dipole moment 


The dipole potential is 


W r ) = 4^ f (?■?') d 3 r' = -^-^v ■ J p(v')v'd 3 r' 


r P 

4neor 2 


(3.7) 


where the dipole moment is 


p = J p( r') r ' d 3 r', 


(3.8) 
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or, for a collection of n point charges located at r^, 

n 

p = 5>rl fe ]. (3.9) 

k =1 

An example is the physical dipole comprising charge +g located at rW and charge — q located 
at rl 2 l for which p = [(+g)rW + (— q)r^] = gd, where d = (rW — rl 2 !) is the separation vector 
from the negative charge to the positive charge. 


3.4 The quadrupole moment tensor 


The quadrupole potential is 


^Wr) = Jp( r ')( r ') 2 \ [ 3 (r • r ') 2 - l] dV. 


(3.10) 


Exercise 1-2 is to show that this can be written 


1 1 


3 3 


r q uad( r ) - 47r£or 5 2 


(3.11) 


i=l j =1 


where rq = x, = y and r% = z, and the quadrupole moment tensor for charge density p( r), 
or for n charges with the fcth charge located at r^, is 


Qij — 


[ p(r') [Sr'r-- Sij (r 1 ) 2 ] d?r', Qij = ^q k 3rf ] rf ] -5ij(r [k] ) 2 . (3.12) 

J k =1 


As an example, consider the physical quadrupole made up of four charges at the corners of a 
square of side 2 m in the xy plane: + 10 -4 C at ( 1 , 1 , 0 ), - 10 -4 C at (- 1 , 1 , 0 ), + 10 -4 C at (- 1 ,- 1 , 0 ) 
and -10 -4 C at (1,-1,0), with distances in metres. The charge arrangement is depicted below 
together with its quadrupole moment tensor 


1 

o 

T —1 

1 

+ 10“ 4 

+1CT 4 

-lcr 4 


Qij — 


0 1.2 x 10“ 3 0 

1.2 x 10" 3 0 0 

0 0 0 


(C m 2 ). 
(3.13) 
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Since the monopole and dipole terms are zero, the quadrupole moment tensor is independent of 
choice of origin. The potential calculated using the quadrupole moment is plotted in Fig. 3.2, 
where it is compared with the potential calculated for the four charges directly using Coulomb’s 
law. Another example of the quadrupole moment tensor is given in Exercise 3-4. 



Figure 3.2: Potential due to the four charges and their quadrupole moment tensor (Eq. 3.13) 
as described above: left — calculated using the quadrupole moment; right — calculated for the 
four charges directly from Coulomb’s law. (Potentials higher than 10 6,4 V are set to 10 6,4 V). 
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3.5 Electric dipoles 


Electric dipoles are important because, as we shall see in the next chapter, when dielectrics are 
subject to an electric field they acquire an electric dipole moment. Also, oscillating dipoles are 
efficient radiators of electromagnetic waves, and are discussed in the 2nd book in this series 
“Essential Electrodynamics”. 


For a dipole p at the origin pointing in the ^-direction 


^dip(r) = 


1 p • r 1 p cos 9 


47T£o r 2 47T£o r 2 


(3.14) 


The components of the electric field, Edi p (r) = —VVdi p (r), are 


_ 8V 2 p cos 9 1 8V p sin 9 

Er = -^ = — - E e = —~- = f- 


dr 4:7T6or 3 


r 89 Amor 3 


1 8V 
rsin0 8(f) 


= 0 . 


Edi P (r) = 


P 


47T£o^ 3 


2cos0r + sin 9 6 


(3.15) 


For a dipole at the origin pointing in an arbitrary direction 


Edip(r) = - V 


1 p • r 


47T£o J ’ 


(3.16) 


Edip(r) = 


47T£q 


47T£q 


47T£o^ 3 


P-rV(4V4v(p-r) 


P' r ( -3 ) ( ~ 4 p + ~ 3 V {PxX + p y y + p z z) 


[3(p • f)f — p] ■ 


(3.17) 

(3.18) 

(3.19) 


While a physical dipole comprises two charges, +q and —q, separated by distance d and has 
dipole moment p = qd where d extends from — q to +q, a pure dipole can be thought of as a 
physical dipole with dipole moment p = gd, with q = p/d, in the limit that d —>► 0. Field lines 
of a physical dipole and a pure dipole are shown in Fig. 3.3. If we were to zoom out so that the 
scale was very much larger than the separation between the two point charges in part (a) the 
two fields would then appear to be identical. 
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(a) (b) 




Figure 3.3: Field lines of (a) a physical dipole, and (b) a pure dipole. 

3.6 Multipole expansion in spherical coordinates 


In addition to the multipole expansion in Cartesian coordinates, one can perform a multipole 
expansion of a localised charge distribution in spherical coordinates. We give this without proof, 


V(r,e,<j>) = 52 o7TT r (£+1)l WM), 


£(3 z — J z — J T 1 

U £=0 m——i 


(3.20) 


where the spherical multipole moments are 
qt,m = j Yl m (6', (f)')r ,e p(r')d 3 r'. 

These are directly related to the Cartesian multipole moments, e.g. 

1 / 3 [ 3 

goo = ^<?, qn = -^(Px-ip y ), < ?io = y^P- 

Summary of Chapter 3 


(3.21) 


(3.22) 


Multipole expansion for a localised charge distribution 
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— Far from a localised charge distribution the potential is 


V(r) 


q r • p 11 ^ ^ 

+ --o + --^ / v Qij r i r j + 


Aive^r Attest 2 Attsq r 5 2 


i=l j =1 


(3.23) 


— Dipole moment for charge density p{ r), and point charges 


~ n 

p= I p(r')r'd 3 r', p = ^C®r [?] (3.24) 

2=1 

Quadrupole moment tensor for charge density p( r) and point charges, 


Q 


13 


= [ p( r ') [3 r{rj - <% (r') 2 ] dV, Q {j = ^q k 3rf ] rf ] - <% (r^) 2 (3.25) 

J k = 1 


Potential and electric field of a dipole p at the origin pointing in the ^-direction 


= Edip(r) = P 




2 cos 6 r + sin 9 0 


(3.26) 
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Exercises on Chapter 3 


3-1 On the surface of a non-conducting sphere of radius a is surface charge density cr(a, 0, 0) 
ao cos 3 0. Find the dipole moment of the sphere. 


3-2 The quadrupole potential is 



J [3(r-r ') 2 - 1 ] d 3 r' 


(3.27) 


Show that it can be written as 



i=l j =1 


(3.28) 


where the quadrupole moment tensor is 


Qij 


J P ( r 0 [3 



% (r 1 ) 2 ] d 3 r'. 


(3.29) 


[This exercise is easy using index notation.] 

3-3 A physical quadrupole is made up of four charges lined up along the 2 axis: -qo at ( 0 , 0 ,— 2 a), 
+go at (0,0,—a), +go at (0,0,a) and -qo at (0,0,2a). (a) Obtain the quadrupole moment, 
(b) Find the potential at r = ( 6 , 6 , 0) for b >> a. 

3-4 Charge — q is located at the origin and charge +q is located at 

(x,y,z) — (a sin 0 o cos 0 o, a sin 0 o sin 0 o, a cos 0 o)- 

(a) Find the the non-zero moments of the multipole expansion of the potential in Cartesian 
coordinates, i.e. < 7 , p, Qij (if non-zero), and use these moments in the multipole expansion 
in Cartesian coordinates to find the potential at (x, y, z) = (r sin 0 cos 0, r sin 0 sin 0, r cos 0) 
where r a. (b) Find the non-zero moments of the multipole expansion of the potential 
in spherical coordinates, i.e. 



(3.30) 


and use these moments in the multipole expansion in spherical coordinates to find the 


potential at (r, 0, 0) where r a. Compare the result with that from part (a). 
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4 Macroscopic and microscopic dielectric theory 

Learning objectives 

— To understand that dielectric materials when placed in an electric field acquire an electric 
dipole moment per unit volume described by the polarisation field, and that this leads to 
polarisation volume and surface charge densities. 

— To learn how Gauss’ law must be modified to include polarisation charges, and that it 
may be written in terms of a new field called the electric displacement. 

— To learn that for many dielectrics the polarisation is proportional to the electric field 
with the constant of proportionality being the electric susceptibility, and to learn how the 
susceptibility is related to the permittivity. 

— To be able to calculate the capacitance of capacitors having a simple geometry. 

— To know and be able to apply the boundary conditions on the potential, electric field and 
displacement field at an interface between two dielectrics. 

— To understand, on a microscopic level, the causes of polarisation, including the orien¬ 
tational polarisability of polar molecules, electronic polarisability of non-polar molecules 
and the ionic polarisability of crystals. 

— To know that the electric field at the atomic scale responsible for aligning or stretching 
atomic or molecular dipoles is different to the macroscopic electric field in the dielec¬ 
tric, and that the Clausius-Mossotti formula accounts for this in relating the molecular 
polarisability to the permittivity. 


4.1 Dielectrics 


Electrons in dielectrics or insulators are attached to specific atoms or molecules. This is very 
different to in conductors where some fraction of electrons are free to move throughout the 
material. Some molecules have permanent electric dipoles. If an electric field is applied to a 
dielectric, molecules rotate or atoms become stretched, and the material becomes polarised. A 
polarised material has a net dipole moment per unit volume called the polarisation field , P(r) 
(units: C m -2 ). 
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4.2 Macroscopic dielectric theory 


By working out the potential of a polarised object we shall identify parts due to volume and 
surface polarisation charges. The electrostatic potential due to a point dipole located at r' is 


W r ) 


1 p • R 

4 TTSo B? 


(4.1) 


where R = (r — r') as usual. For a polarised dielectric the dipole moment due to volume element 
d 3 r' at position r' is P(r / )d 3 r / , and so for the whole object of volume V 


nr) 


J_ f *V) • ^ 

4vt£ 0 ,/v i? 2 



(4.2) 


But we know (Eq. B.9) that R 2 R = V'i? 1 where V 7 = 


\ 8 ^ d ^ d \ 

x & ; + y 57 +z d?J 


and so 


V(r) = 



(4.3) 
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Then we can use the product rule V ' • (uA) — (V r u) • A + u( V' • A) to integrate by parts 


E(r) = —— [ V'. [ 

Attsq J v R 47 T£ 0 J v 


r , 

R dr ' 


(4.4) 


and then use Gauss’ theorem on the 1st integral to obtain 


F(r) = 


1 _ 1 /V' P(O j3 „, 


4ttsq 


R 


4ttsq 


L 


R 


-dV, 


(4.5) 


where surface S bounds volume V. 

Comparing Eq. 4.5 with the potential due to surface and volume charge densities (Eqs. 1.20) 
we see that 


V(r) = -—<f ^dS' +-*— I ^^d 3 r\ 
4tt£o Js R 4tT£q J v R 


(4.6) 


and the “polarisation” or “bound” volume and surface charge densities are 


Ppoi(r) = — V • P(r) (unit: Cm 3 ), 


crpoi(r) = P(r) • n (unit: C 


m 


- 2 ^ 


(4.7) 

(4.8) 


Figure 4.1 is a cartoon which attempts to illustrate how surface and volume polarisation charges 
occur at an atomic level in a dielectric material. In Fig. 4.1(a) we show a microscopic rectangular 
slab of dielectric of length L containing an array of atoms which have been stretched to form 
tiny dipoles in the presence of an applied electric field. The material is uniformly polarised, 
with P(r) = Pox. On the right hand side of the dielectric cr po i(L, y,z) = P • x = +Pq (positive 
surface charge), whereas on the left hand side cr po i(0, y, z) — P • (—x) = —Po (negative surface 
charge). There is no surface polarisation charge present on the top and bottom surfaces as 
P • n = 0 there. Also, because P is uniform, p po i(r) = —V • P(r) =0 and so there is no volume 
polarisation charge either. 

In Fig. 4.1(b) the dielectric has been polarised non-uniformly and in the rectangular slab shown 
P(r) = (Po + ax) Sc. There is still negative surface charge on the left hand side and positive 
surface charge on the right hand side, but they have unequal magnitudes: cr po i(0, ?/, z) = —Po 
and cr po i(L, y, z) = (Po + aL). Inside the slab V • P(r) = a and so the polarisation volume 
charge density is negative, p po \(r) = —a. This can be expected as, although there are equal 
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(a) P(r) is uniform (b) P(r) increases with increasing x 



Figure 4.1: Explanation of polarisation charge: (a) uniform polarisation giving rise to surface 
polarisation charge on the surfaces defined by x = 0 and x = L, (b) non-uniform polarisa¬ 
tion giving rise also to a volume polarisation charge density. The displacement of electron 
distributions from the nuclei are exaggerated. 


numbers of real positive and negative charges, the negative charges occupy a smaller volume 
than the positive charges (compare the two boxes in Fig. 4.1b) giving rise to a net volume 
polarisation charge density which is negative. If the slab has cross-sectional area S, then the 
total polarisation charge is q po \ = S[(—Po) + (Po + °P) + (—a x L)\ = 0. 


4.3 Gauss' law and the electric displacement 


Inside dielectrics, Gauss’ Law must include both free charge pf and polarisation charge p po \ 


, E ( r ) = M r ) + i P oi( r ) = M r ) 

£0 


V • P(r) 


£o 


V-[e 0 E(r)+P(r)] = p f ( r) 


(4.9) 

(4.10) 


Hence, Gauss’ law in matter becomes 

V-D(r) = p f { r), fn-dS = Q f , (4.11) 

where D(r) = [£oE(r) + P(r)] is the displacement field (unit: C m -2 ). 
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4.3.1 Energy density 


We shall now find the energy density of an electric field in a dielectric, starting with 


W 


= \ f pf{r)V(r)d 3 r, 

^ J all space 


and using Gaus’ law to replace pf 


W = - 


\\ (V-D )Vd 3 r, 

^ Jail space 

= 11 V • (VO)d 3 r - \ [ (W) • D d 3 r, 

^ J all space ^ J all space 


(4.12) 


(4.13) 

(4.14) 


where we have used a product rule to integrate by parts. Then using Gauss’ theorem on the 
1st integral and E = — 'VV on the 2nd integral we find 



S'at oo 


(l/D) • dS + 


\L 


E • D d 3 r, 


all space 


(4.15) 
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all space 


E • D d 3 r, 


(4.16) 


since the surface integral is zero because for a localised charge distribution VD ~ r 3 and 
S ~ r 2 as r oo. Hence the energy density is ^E • D. 


4.3.2 Linear dielectrics 


For linear dielectrics, i.e. most materials unless the electric field is too high, 


P = £0XeE 


(4.17) 


where Xe is the electric susceptibility. Then 


D — (<sqE T P) — £q(1 T Xe)E — <sE — s r s qE — Ks qE 


(4.18) 


e is the permittivity and e r — K is its relative permittivity or dielectric constant. Some values 
of relative permittivity are: e r — 1 (vacuum), e r =l .00055 (Nitrogen), £ r =5.9 (salt), £ r =80.1 
(water). 

4.3.3 Capacitors 

A capacitor is a device for storing charge, typically comprising two conductors separated from 
each other by air or a dielectric, the simplest being the parallel plate capacitor comprising 
two parallel conducting plates separated by a small gap. The potential difference between two 
conductors, “a” and “b”, say, is 



(4.19) 


Jb 

The electric field due to the two charged conductors which have surface areas S a and Ss is 



(4.20) 
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where <r(r) is the surface charge density on the conductors. If the charge on each conductor is 
doubled, E is doubled and V must also be doubled, and so the charge stored q is proportional 
to V. The constant of proportionality is the capacitance C, such that q = CV. The unit of 
capacitance is the farad, and is named after Michael Faraday. 


To charge a capacitor, you must transport electrons from the positive plate to the negative 
plate. If the capacitor has charge q ', the potential difference is V = q'/C, and so the work done 
to add a small additional charge dq' to a capacitor which already stores charge q' is 


dw = V dq 1 = ^ dq (4.21) 

We can integrate this from q' = 0 to q' = q to determine the work done to charge an initially 
uncharged capacitor with charge q 


W = 





(4.22) 


A schematic diagram of a parallel plate capacitor is shown in Fig. 4.2(a) where the plate 
separation is greatly exaggerated. The capacitance can be calculated as follows: (i) neglecting 
fringing fields (bulging out around the edges of the plates) the electric field is located solely 
between the plates and is perpendicular to the plates, (ii) using Gauss’ law the surface charge 
density is ao = £qE, (iii) the total charge is 

q = Sa 0 = Se 0 E = Se 0 — . C = —. (4.23) 

a a 


If a dielectric is placed between the plates, E is unchanged (it depends only on V and a), but the 
dielectric becomes polarised P = EoXeE and acquires surface polarisation charge (see Fig. 4.2b) 

(Jpoi = p n = Xe^oE • n. (4.24) 

To keep E the same, more charge (shown in red) must flow onto the plates to cancel the 
polarisation charge such that 

cr = (cro + crp 0 i) = (sqE + 6ox e E) = (l + Xe)^o^ = zE. .’. C — (4.25) 
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° 0 =£ 0 E 




Figure 4.2: Parallel plate capacitor filled with: (a) air (e « so), (b) dielectric with permittivity 
e. The separation between the plates is greatly exaggerated, and the gap between the dielectric 
and the plates is negligible. 


Practical capacitors used in ordinary circuits are made of two long (to make S large) strips of 
metal foil separated by a thin (to make a small) insulating layer, and rolled up into a cylinder 
(to save space). 
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4.4 Boundary conditions on V, E and D 


The boundary condition on V at the interface between two materials is simply that V must 
be continuous across the interface; otherwise at the interface the electric field would be infinite 
there. 

To obtain the boundary conditions on E and D at the interface between two materials we split 
up the problem into finding the conditions on the components normal to and parallel to the 
interface. We obtain the boundary conditions on the normal components by applying Gauss’ 
law to an infinitessimally-thin Gaussian pill-box spanning the interface between two materials 
as in Fig. 4.3(a). Since both the surface polarisation charge and the free surface charge will 
enter into Gauss’ law for E, whereas only the free surface charge enters into Gauss’ law for D, 
the boundary condition on the normal component of D is much more useful. Then, 


(p D(r) • n dS — Qf , e neb Di • ni£i + D2 • 112S2 = Si a /. 

JS1+S2 


(4.26) 


(Di - D 2 ) • ni = a f , 


(4.27) 


since S 2 = £ 1 , an d so the normal component of D is unchanged across the interface between 
two materials, unless there is free surface charge on the boundary. 



Figure 4.3: (a) Gaussian pill box used to find the boundary condition on the normal component 
of D; (b) integration loops used to find the boundary condition on the parallel component of 

E. 


In electrostatics the electric field is conservative, and so 


1 


E(r) • dr = 0 


(4.28) 


for any “loop”. Using infinitesimally narrow loops T a or T^ of lengths £ a and <4 spanning the 
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interface, and with t a orthogonal to t&, as in Fig. 4.3(b), we obtain 

(Ei - E 2 ) • (4t„) = 0, and (E, - E 2 ) • (44) = 0. (4.29) 

Adding the components of E parallel to t a and to t& we have the component of the electric field 
parallel to the interface, Ell, and this must be the same on either side of the boundary, 

e! = e|. (4.30) 

4.5 Microscopic dielectric theory 


In the presence of an electric field, atoms become stretched or molecules rotate thereby acquiring 
on average a dipole moment p = <a po iE where <a po i is the polarisability. 


4.5.1 Orientational polarisability of polar molecules 

Molecules having a permanent electric dipole are called polar molecules. A common example 
is the water molecule illustrated in Fig. 4.4(a), but there are many other examples including 
carbon monoxide and ozone. In the presence of a uniform electric field the net force on an 
electric dipole is zero, but a permanent electric dipole such as a water molecule will experience 
a torque tending to align its dipole moment with the electric field (Fig. 4.4b), 

N = (F) x (gE) + (-|) x (-gE) = gL x E, 

N = p x E. (4.31) 

At a finite temperature, collisions with other molecules regularly mis-align individual molecules 
such that, on average, only partial alignment is maintained. 

The lowest energy state of a polar molecule in an electric field corresponds to when p is aligned 
parallel to E, and the work done to rotate p so that it is oriented at angle 6 to the electric field 
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Figure 4.4: (a) Diagram of water molecule showing its dipole moment, (b) Forces on individual 
charges of a physical dipole. 


is 


0 rO 

N(6')d9' = / pEsm6'd6' = pE( 1-cos#). (4.32) 

Jo 

At finite temperatures, the probability density of W, and consequently of cos#, can be obtained 
using the Boltzmann distribution, i.e. 


W(9) = [ 
Jo 
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f(W) oc e ~ w ^ kT oc e P Ecosd / kT , /(cos#) oc e y 


cos 6 


(4.33) 


where y = pE/kT. For this probability density, the time-average of a dipole moment’s compo¬ 
nent parallel to E is 


(p cos 6) 


f\ cos9 e ycos0 dcos6 

p -;- 

f_ ± ey c °s y dcos6 


= P 


coth(y) 


= pL(y), 


(4.34) 


where L(y) is the “Langevin function” which is plotted in Fig. 4.5. For ordinary electric fields 
and temperatures, y <C 1 such that L(y) ~ y/ 3, and the orientational polarisability obeys the 
“Curie law” named after Pierre Curie (1859-1906) who discovered the 1/T dependence, 


(P-E) „ 

E ~ 3 kT' 


(4.35) 



Figure 4.5: The Langevin function (Eq. 4.34). 


4.5.2 Electronic polarisability of non-polar molecules 


In an attempt to understand how atoms acquire an electric dipole moment when placed in 
an electric field we approximate an atom by a nucleus screened by its inner electrons (effective 
charge +e) and a spherical valence electron “cloud” (charge -e) of radius ao ~ 10 -10 m (Fig. 4.6). 
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In the presence of an applied electric field, as the electron cloud becomes displaced with respect 
to the nucleus by an amount d e it experiences a restoring force 


(d e /q 0 ) 3 e 2 - _ e 

47T£o^e 6 47 T£()Uq 


(4.36) 


This is equivalent to Hooke’s law with “spring constant” k ~ e 2 /47r£o«o- 


No applied field Applied electric field 



Figure 4.6: Semi-classical model of an atom without, and with, an applied electric field (dis¬ 
placement greatly exaggerated). 


The net force on the valence electron is then 


m e 


h 2 de 

dt 2 


—eE — ftd e , 


(4.37) 


which must be zero in electrostatics. Hence there will be an equilibrium displacement, dipole 
moment and molecular polarisability 


E , 2 E 

= ~e— , p = -ed e = e —, 


_ P 

a P° 1 — £ 


e 

K 


= 47T£ 0 «0- 


(4.38) 


4.5.3 Ionic polarisability of crystals 

Consider a tiny piece of a salt crystal sketched in Fig. 4.7(a) comprising Na + and Cl - ions in 
a face-centred cubic structure. In the absence of an applied electric field the crystal has zero 
net dipole moment. Within the crystal, physical dipoles are associated with pairs of Na + and 
Cl - ions but they cancel out over the volume of the crystal as shown (there are equal numbers 
of dipoles pointing to the left as there are to the right). This might look a bit suspicious, as 
you could draw all the dipoles to point to the right as in Fig. 4.7(b). However, for this rather 
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artificial construct there will remain unpaired ions, resulting in negative surface charge on the 
right and positive surface charge on the left, and the dipoles associated with the surface charges 
would cancel the internal dipoles. 



564 pm 


Figure 4.7: Face-centred cubic crystal structure of NaCI (salt) showing: (a) locations of Na + 
and Cl - ions and how dipoles associated with ion pairs will cancel over the volume, (b) an 
alternative way of showing the cancellation as in this case internal dipoles are cancelled by 
those associated with surface charge, (c) a net dipole moment if an electric field is applied (ion 
displacements exaggerated). 


Fig. 4.7(c) shows negative ions displaced to the left and positive ions displaced to the right in 
response to an applied electric field such that the internal dipoles no longer cancel, but give 
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rise to an overall polarisation in the direction of the applied field. Usually for crystals (NaCl is 
an exception) the direction of E relative to crystal planes matters, and p is related to E by the 
polarisability tensor, rather than the (scaler) polarisability. 


4.5.4 Clausius-Mossotti formula 

In an amorphous substance made of non-polar molecules, or an isotropic crystal such as NaCl, 
placed in an external electric field E ext a molecule located at r mo i will acquire an induced dipole 
moment 


P — ^polPmicro (^mol) 


(4.39) 


where E m i cro (r mo i) is the microscopic electric field at the location of the molecule. 

In a slab of dielectric in a uniform external electric field E ext , perpendicular to the slab, the 
macroscopic field inside the slab will be 

E ma cro = —E ex t- (4.40) 

£ 

But this is not the electric field we need for Eq. 4.39 — we need the field that would be at r mo i 
if the molecule weren’t there. 

For amorphous substances we can approximate the microscopic electric field by the field at the 
centre of a spherical hole of radius a and volume ~ 1/N surrounding the molecule, where N is 
the number density of molecules (Fig. 4.8). At the centre of the hole, u O”, we need to add the 
field due to the surface polarisation charge of density cr po i(0) on the surface S of the hole, 

E mtao = E,„„ ro + -j-/ = E macro + -i- / dS. (4.41) 

4tt£o Is R 4tt£o J s a 2 

The surface polarisation charge density is 

<7poi(0) = P-n = P z-n —PcosO. (4.42) 

The minus sign comes from the direction of n which is outward from the volume of the material 
(towards centre of hole). 
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Figure 4.8: Geometry used for calculating microscopic electric field inside an amorphous di¬ 
electric. 


From symmetry arguments E m i cro must be in the z direction, and using n • z = — cos 6 we have 


E m icro ‘ z — F/ macro + 


1 


Eri 


47T£:oa 2 J_i 

P 


(—P cos 9)( — cos 9)27ra 2 d(cos 6 ), 


Z F^macro + “— / COS 6 d COS 6. 


/: 


E m icro — E macro 


+ 


2eo J -i 

P 

3e 0 ’ 


(4.43) 

(4.44) 

(4.45) 


We are now in a position to obtain the Clausius-Mossotti formula, named after German physicist 
Rudolf Julius Emanuel Clausius (1822-1888) and Italian physicist Ottaviano-Fabrizio Mossotti 
(1791-1863), which relates the molecular polarisability to the relative permittivity. If the num¬ 
ber density of molecules is N, the polarisation will be 


P — Na pol E m j Cro 


Na po \ I E 



(4.46) 


Solving for P, 


P = 


Not po l 


1 Not poi/3^0 


E„ 


Hence, we obtain the susceptibility and the relative permittivity 


Xe = 


P Na pol / £q 

e 0 E 1 - lVa po i/3£ 0 ’ 

1 + 2A^a po i/3£o 


£ r — (1 + Xe) — 


1 - Na p o i/3e 0 ’ 


(4.47) 


(4.48) 

(4.49) 
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and finally the Clausius-Mossotti formula 


C^pol 


3^0 / s r - 1 \ 

N \e r + 2J' 


An example is given in Exercise 4-7 


(4.50) 


Summary of Chapter 4 

Macroscopic dielectric theory 

— The polarisation field P (units: C m -2 ) is the net dipole moment per unit volume. 

— Surface and volume polarisation charge due to polarisation of dielectric 

Ppoi(r) = —V • P(r) (unit: C m -3 ), cr po i(r) = P(r) • n (unit: C m -2 ). (4.51) 
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— Gauss’ law and the electric displacement field D(r) = [eroE(r) + P(r) 


V • D(r) = Pf(r), D • dS = (5/, enc l- 


(4.52) 


— For linear dielectrics P = £oX e E 

D = (^oE + P) = 6?o(l Xe)E = £~E = £ r £gE = (4.53) 

where e is the permittivity, K — e r is the relative permittivity or dielectric constant. 

— Boundary conditions at interface between two dielectrics 

(Di — D 2 ) • ni = <Tf, E| = e||. (4.54) 


Microscopic dielectric theory 


— In an electric field, atoms or molecules acquire a dipole moment p = oE m i cro where <a po i 
is the molecular polarisability. 

— The “microscopic electric field” at the location of an atom or molecule differs from the 
macroscopic field in a dielectric. In an amorphous dielectric 

P 

E m icro = E macro + — 

— Clausius-Mossotti formula 

3^o ( e r - 1 \ 

“ pol_ N (e r + 2 )■ 

Exercises on Chapter 4 

4-1 A dielectric sphere (dielectric constant K ) of radius a is placed in an initially uniform 
electric field Eo- (a) What are the boundary conditions on V, E and D for this problem, 
(b) Find the potential everywhere, (c) Find E, D and P everywhere, (d) Find the dipole 
moment of the sphere and the surface polarisation charge density. 


(4.55) 


(4.56) 
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4-2 An electret, i.e. a piece of material with a permanent electric polarisation, is in the shape 
of a sphere of radius a and has P(r) = Po- (a) Find the surface polarisation charge density 
and the dipole moment of the sphere, (b) find V, E, and D everywhere, and (c) sketch the 
field lines of E and D. 

4-3 The space between two concentric conducting cylinders of radius a and b > a and length 
L b is filled with a dielectric with permittivity e. The inner and outer conductors 
are held at potentials V a and V&, respectively. Find: (a) E, D and P everywhere; (b) 
the polarisation surface and volume charge density everywhere, and the net polarisation 
charge; (c) the free charge on the inner and outer conductors, and the capacitance. 

4-4 A spherical capacitor is filled with two different dielectrics with permittivities s\ and e 2 as 
shown in the diagram. The capacitor is charged such that charge -\-q is on the inner con¬ 
ductor. Find: (a) D, E and P everywhere; (b) the polarisation surface and volume charge 
density everywhere; (c) the net polarisation charge; (d) the potential difference between 
the inner and outer conductor, and the capacitance of the capacitor. 



4-5 A uniform slab of material with relative permittivity e r is suspended parallel to the xy- 
plane, and has its lower surface at z = 0 and its upper surface at z — d. Outside the slab 
there is a uniform electric field Eg = Eq (sin 6q x — cos 9q z) . (a) Find formulae for E, D and 
P in the dielectric, the angle between E in the dielectric and the normal to the surface, 
and the surface polarisation charge density at 2 : = 0 and z — d. (b) Find numerical values 
for the case of Eo = 1000 V m -1 , 6q — 45° and e r — 2, and include a sketch showing field 
directions. 

4-6 Derive the force on an electric dipole in a non-uniform electric field. 

4-7 The relative permittivities of Nitrogen, Argon and Hydrogen in gas (at 20°C) and liquid 
phases are given below. 
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Element 

N 2 gas 

Ar gas 

H 2 gas 

N 2 liquid 

Ar liquid 

H 2 liquid 

s r 

1.000546 

1.000517 

1.000272 

1.45 

1.53 

1.22 


http: / / www.kayelaby.npl.co.uk/general_physics / 22 / 2_2_1 .html 


Use the Clausius-Mossotti formula to find the electronic polarisability, and compare the 
results for the same elements in the liquid and gas phases. [You will need to look up any 
constants and the atomic weights and densities required.] 
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5 Magnetic field and vector potential 

Learning objectives 

— To understand the concept of the magnetic field, and that electric currents cause magnetic 
fields. 

— To learn that an electric charge moving in a magnetic field experiences a force given by 
the Lorentz force law, and that as a result wires carrying electric currents in a magnetic 
field experience a force given by Ampere’s force law. 

— To be able to apply the Biot-Savart law to calculate magnetic fields due to simple current 
distributions. 

— To be able to calculate the vector potential for simple current distributions, and to use it 
to find the magnetic field. 

— To be able to use Ampere’s law in integral and differential form to find the current density 
from magnetic fields, and the magnetic fields from symmetric current distributions. 

— To understand the concept of the magnetic flux, that it must be zero for a closed surface 
(no magnetic charge law), and how the magnetic flux linking two circuits leads to the 
concept of mutual inductance. 

— To learn about the multipole expansion of the vector potential and that the dipole term 
is the lowest non-zero term. 

— To be able to obtain the magnetic dipole moment of a current distribution, and the 
magnetic field of a magnetic dipole. 

— To be able to calculate the torque on a magnetic dipole or a circuit in a magnetic field. 


5.1 Magnetism 


Magnets have been known since antiquity - pieces of lodestone (magnetised magnetite FeaC^) 
attract pieces of iron, and when suspended freely point to the magnetic North Pole. Magnetite is 
a ferromagnetic mineral which became magnetised by strong magnetic fields, possibly associated 
with lightning strike. Compasses were invented in China in the 2nd Century B.C. and were 
being used by European navigators by the 12th Century A.D. A compass points in the direction 
of the magnetic field. 
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Magnetic fields occur throughout the Universe - the Earth’s magnetic field is probably produced 
by electrical currents associated with magnetohydrodynamic flows within the Earth’s core. 
Similarly, the Sun and the stars have magnetic fields, pulsars being extreme examples, with 
some varieties known as magnetars having fields as high as 10 14 times that on the Earth’s 
surface. Magnetic fields also pervade the Solar System, interstellar space and, at some level, 
everywhere in the Universe with the possible exception of giant voids in cosmic structure. 

Apart from permanent magnets, the main sources of magnetic fields are electrical currents. In 
the case of a current along a long straight wire, the magnetic field circulates around the current 
as shown in Fig. 5.1(a). The unit of magnetic field, B, is the tesla (T) which is named after 
Serbian-American electrical engineer Nicola Tesla (1856-1943). 


5.2 The magnetic field and forces 


As well as producing magnetic fields, electrical currents experience a force in the presence of 
a magnetic field. For example, the current in one wire produces a magnetic field which causes 
the current in another wire to experience a force. Historically, this magnetic force between two 
wires carrying currents was investigated in the 1820’s and 1830’s by Gauss and French physicist 
Andre-Marie Ampere (1775-1836). A formulation giving the force on a particle with charge Q 
moving with velocity v through an electric field E and magnetic field B is the Lorentz force 
(see Fig. 5.1b), 

F = Q(E + v x B), (5.1) 

named after Dutch physicist Hendrik Antoon Lorentz (1853-1928), of Lorentz transformation 
fame. 


5.2.1 Forces on currents 

We can consider the force between two parallel wires using the Lorentz force law as in Fig. 5.1(c). 
We define a line current I as the net charge per unit time flowing along a filamentary wire, and 
so the net charge flowing past some point in time dt is dQ — Idt. The unit of current is the 
amp (A) and is named after Ampere. Note that 1 A = 1 C s -1 . 

If the effective current flow speed is v then charge dQ is located over length dr = v dt of wire, 
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Figure 5.1: (a) Circulation of magnetic field around a long straight current, (b) Lorentz force 
on charge Q moving at velocity v parallel to the wire carrying current I. (c) Lorentz force on 
length dr of a second wire (carrying current I 2 ) due to the magnetic field of current I\. 


and hence the force on section dr is 


rfF mag = dQv x B = (Idt)w x B = Idr x B. 


(5.2) 


Since the force on section dr of the wire is dF mag = I dr /. B the force on the entire circuit is 
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F mag = (Idr x B). (5.3) 

This is one version of Ampere’s force law. 

Currents can also flow across a surface as “current sheets” and are described by the surface 
current density K (amps per metre). Currents flowing through a volume are described by the 
volume current density J (amps per square metre). The latter two cases are illustrated in 
Fig. 5.2. Since the force on section dr of the wire is dF mag = I dr x B the forces on the entire 
circuit for the three cases are: 



Figure 5.2: Sketch illustrating (a) surface current density (b) volume current density. 


5.2.2 Hall effect 

As discovered by American physicist Edwin H. Hall (1855-1938), when a current is passed 
through a semiconductor in the presence of a magnetic field a voltage, called the Hall voltage , 
develops perpendicular to the current direction. This phenomena, due to the Lorentz force on 
electrons or holes carrying the current, is called the Hall effect. The geometry for calculating 
the Hall voltage is shown in Fig. 5.3. 

Note that the direction of the current / in the diagram is that of the conventional current, so 
that the motion of electrons is in the opposite direction. The current is 

/ = Jhd — (n c qv) hd, (5.5) 

where J is the current density, n c the number density of charge carriers of change g, and v their 
effective velocity. The Lorentz force gvx B acts on the charge carriers deflecting them towards 
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Figure 5.3: Geometry for calculating the Hall voltage. 


one side, and causing one side of the semiconducting strip to become positive and the other 
negative. 


In equilibrium, the Lorentz force of magnitude qvB is opposed by the electrostatic force of 
magnitude qVn/d due to this charge separation, where Vh is the Hall voltage. Hence, 


Vh — dvB — d 


I 1 IB „ IB 

- = - T = Rh ~T' 

n c qhd n c q h h 


(5.6) 


The ratio of the Hall voltage to the product of the current and the magnetic field divided by 
the element thickness is known as the Hall coefficient, 


Rh = 


n c q 


(5.7) 


It is a characteristic of the material from which the semiconductor is made, as its value depends 
on the type, number and properties of the charge carriers: electrons in a conductor (q — — e), 
electrons in a “n-type semiconductor” (q = —e), and holes which act like positive electrons in 
a “p-type semiconductor” (q = +e). 


The Hall effect can be used to measure magnetic fields with a “Hall probe”. The quantum 
Hall effect which occurs for a two-dimensional electron gas in high magnetic fields and low 
temperatures was discovered by in 1980 by Klaus von Klitzing (b. 1943) and has led to a new 
standard for measurement of electrical resistance. 


5.2.3 Cyclotron motion 

A charged particle with velocity perpendicular to a uniform magnetic field will move in a circular 
orbit as in Fig. 5.4. The radius of the orbit and its angular frequency can be found by equating 
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the Lorentz force QvB to the centripetal force mv 2 /R, giving for non-relativistic velocities 


R = 


mv 
QB ’ 


v _QB 

^cycl j 

R m 


(5.8) 


where cu cyc \ is called the cyclotron frequency. For relativistic particles of total energy 7771 c 2 and 
momentum 7 mv 



Figure 5.4: Circular motion of a charged particle in a magnetic field. 

If v is not perpendicular to B, the particle will undergo helical motion. 
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5.2.4 Biot-Savart law 

Named after French physicists Jean-Baptiste Biot (1774-1862) and Felix Savart (1791-1841) 
this experiment-based law gives the magnetic field due to arbitrary line currents (Fig. 5.5a), 
surface currents and volume current distributions (Fig. 5.5b) 


B(r) = 


po I f dr' x R 


47r 


R 2 


<—> 


Mo 

47r 


K(r') x R 

w 


dS' > 


Mo 

47T 


/ 


J(r') x R j3 , 


R 2 


dr, 


(5.10) 


where R = (r — r') and /j,q — An x 10 7 H m 1 is the permeability of free space — the henry 
(H) is the unit of inductance and is named after American scientist Joseph Henry (1797-1878). 

(b) 




Figure 5.5: Contribution dB of (a) section dr' of line current I, and (b) volume current density 
J(r') inside volume d 3 r', to magnetic field at r. 


We can use the Biot-Savart law to find the magnetic field of a long straight wire carrying current 
/ as in Fig. 5.6. From the cross product, MB must be into the page as shown, i.e. in the <p 
direction in cylindrical coordinates (p, The contribution to the magnetic field at (p, </>, 0) 
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due to the section of wire of length dz' at z' is 


0 ) = 


po I dz' x R 
47T W ' 


p>/ , n x Mo-f f 00 sin(9 r - 




Pol 

4n 

Pol 

47r '"y.oo (p 2 + z /2 ) 3 / 2 


— OO 

’OO 


/ 


dz 7 </>, 


B(p,0,O) = ^p 
Pol 


+ Z ‘ 


/2 


(f). 


• B(p,0,z) = 


27rp 


0 


(5.11) 

(5.12) 

(5.13) 

(5.14) 

(5.15) 

(5.16) 


since the magnitude can only depend on p and I, and B is in the 4> direction as in Fig. 5.1(a). 



Figure 5.6: Geometry for calculation of the magnetic field due to a long straight wire carrying 
current /. 


5.3 The magnetic vector potential 


From Biot and Savart’s law we can find a potential which is a vector field, called the magnetic 
vector potential A(r), from which we can obtain the magnetic field B(r) in an analogous (but 
different) way to obtaining the electric field from the electrostatic potential. We start with the 
Biot-Savart law, and use an identity for Vi? -1 (Eq. B.7), 


B(r) = 


Po 

47T 


J(r') x 


R 


d 3 r' = 



(5.17) 
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Then, using the product rule V x (uF) = (Vu) x F + u(V x F) we have 


"J(r') x V ( — 


= V x 


J(Q 

R 


- ^VxJ(r') 


(5.18) 


But V x J(r ') = 0 because here the differentiation is with respect to x. y. z (not x ', y', z'), and 
so 


B(r) = 



(5.19) 


B(r) = V x A(r), 


(5.20) 


where 



is the magnetic vector potential, or just the “vector potential”. 


(5.21) 
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Exercise 5-2 is to use Eq. 5.21 to find the vector potential of an infinite straight wire along the 
z axis carrying current /, although we can actually find this in a simpler way from the magnetic 
field which we know to be /2'Kpcj) (Eq. 5.16), and we will do so now. 


Using B = V x A in cylindrical coordinates 


Mo I ^1_( - f dA p 

2ivp \p dcf) dz ) P \ dz 

Mo^ 9A Z 

2iip dp 

/“* = -£/ 7 

A(r) = -2^1 In (p) + zC 





(5.22) 

(5.23) 

(5.24) 

(5.25) 


where C is an integration constant whose value has no effect on the magnetic field. 


5.3.1 Gauge transformations 

The magnetic field is unchanged by adding the gradient of an arbitrary scalar field U( r) to A, 

A(r) = [ ^Ad z r' + VU{ r), (5.26) 

47r J K 

because V x VC/(r) = 0. The transformation A A + V!7 is called a gauge transformation. 
Exercise 5-3 is to show that 

V • A(r) = W 2 U(r). (5.27) 

The “gauge” is determined by the value of V • A, and in magnetostatics, it is convenient to 
choose V • A = 0 (Coulomb gauge). 
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5.4 Ampere's law 


This law, named after Andre-Marie Ampere, is obtained by taking the curl of B and using the 
2nd derivative rule for curl of curl, followed by V • A = 0 


V x B(r) = V x [V x A(r)] 

= V [V • A(r)] - V 2 A(r 

>o f J(r 7 ' 


= 0- V" 


/ 


An R 


d 3 r' 


Mo 

47T 


/ J (r') v2 (i 


3 ^/ 


d 6 r 


(5.28) 

(5.29) 

(5.30) 

(5.31) 


Note that in the last line the Laplacian operator has moved to the right of J(r') because the 
differentiation here is with respect to unprimed coordinates. Now, Eq. C.39 gives 


V 2 ( 1 ) = -47t 5 3 (R) = -47r5 3 (r-r'). 


V x B(r) = J J(r') [-47r5 3 (r - r')] dV. 

V x B(r) = MoJ(r), 


(5.32) 

(5.33) 

(5.34) 


which is Ampere’s law in differential form. We can use Stokes’ theorem (Appendix C.6) to find 
Ampere’s law also in integral form 


^ B(r) • dr = /xq J J(r) • dS 


(5.35) 


where loop T bounds surface S. 


Ampere’s law in integral form can be exploited in a similar way to Gauss’ law in integral form if 
the symmetry of problem allows. An “Amperian loop” plays an equivalent role to the “Gaussian 
surface” when using Gauss law. As an example, consider the magnetic field of a long straight 
wire carrying current I. There is cylindrical symmetry, and we can for convenience define the 
z axis to be along the wire. Then the magnitude of the field must depend only on cylindrical 
coordinate p and the field must be in the 0 direction as in Fig. 5.1(a). Taking a circular 
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Amperian loop T of radius p as in Fig. 5.7(a) 


<j> B ■ dr = (Jbolench B 2ixp = [4)1, 


(5.36) 


giving 


B= 2^ 


(5.37) 



Figure 5.7: (a) Amperian loop T centred on the axis of a long wire carrying current /. (b) 
A long solenoid (much longer than shown) and a circular Amperian loop Id concentric the 
solenoid’s axis, and a rectangular Amperian loop T 2 with two sides parallel to the solenoid’s 
axis, one of which is inside the solenoid. 
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Another example is the long solenoid as depicted in Fig. 5.7(b). For a circular Amperian loop Ti 
concentric with the axis, I enc \ — I if Ti is outside the solenoid (pi > a), and I enc \ = 0 if it is inside 
(pi < a). Hence, outside the solenoid, the azimuthal component is B ( f ) (p > a,(j),z) — /ioI/2np, 
whereas inside B^p < a,(f),z) — 0. It is widely stated that the magnetic field outside a long 
solenoid is zero, which appears to contradict this result. It is true that the magnetic field outside 
is very much less than inside, and one can design a solenoid with B = 0 outside by having a 
double coil with current flowing out in the -hz direction in one coil, and back in the —z direction 
in the other, such that the net current through the circular loop Ti is zero. 

From the symmetry of the problem, if there were a component of B in the p direction it would be 
the same everywhere over the surface of a concentric cylinder, and this would give § B • dS ^ 0 
which is not allowed if, as it appears, magnetic charge doesn’t exist (as is discussed in the next 
section). We have already found B^ inside and outside the cylinder, so now we examine B z 
using rectangular Amperian loop V 2 with two sides parallel to the solenoid’s axis, one of which 
is inside the solenoid. The current enclosed is / enc 1 = NI, for all P 2 values, so that Ampere’s 
law implies that B Z L , and thus B z , is the same for all p > a. But B z (p — 00 , </>, z) must be zero, 
so that B z (p > a) = 0. Similar arguments show that B z (p < a,(j),z) — £>o, i.e. B z is constant 
inside the solenoid. Ampere’s law for loop V 2 is then 

(f B • dr = polench B zL = ^0^7, (5.38) 


giving 


B (p, (j),z) 


IMynl z inside solenoid 

Hoi / (2np)(j) outside solenoid (usually neglected) 


where n — N/L is the number of coils (or “turns”) per unit length. 


(5.39) 


5.5 Magnetic flux 


The magnetic flux through surface S is defined as 




B(r) • dS 



A(r) • dS, 


(5.40) 
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and using Stokes’ theorem we find that 

A(r) • dr (5.41) 

where loop T bounds surface S. This result is very useful for finding the vector potential from 
the magnetic field in cases where we can exploit the symmetry of the problem, in a similar way 
to finding the magnetic field from the current density using Ampere’s law. 


We shall now use Eq. 5.41 to find the vector potential due to current / flowing in a long solenoid 
(Fig. 5.7b). From Eq. 5.21, one expects there to be a z component of A outside the solenoid due 
to the net current / flowing along the solenoid, but we will neglect this here. Since the current 
in the solenoid circulates in the (f> direction, from symmetry arguments we should expect that 
A will have a (j) component. Applying Eq. 5.41 to loop for a solenoid of radius a as shown, 


27rpi4£(pi,0,z) 


ttpi/iqtiI inside solenoid (pi < a), 
7ra 2 /io nl outside solenoid (pi > a). 


Hence, 


A (p,(/>,z) 


(ponlp/2)<p inside solenoid (p < a), 

(/ion/a 2 /2p)<fi outside solenoid (p > a). 


(5.42) 


(5.43) 


5.5.1 No magnetic charge law 

We can derive this very important law from 

B(r) = V x A(r), (5.44) 

and then using the second-derivative rule V • (V x c) = 0 we obtain, 


V • B = 0. (5.45) 

Comparing this with Gauss’ law in differential form V • E = p/so we see that Eq 5.45 implies 
that magnetic charge does not exist. Some particle theories predict the existence of massive 
magnetic monopoles, but so far none have been found. Magnetic dipoles do exist however, and 
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are very important in discussing magnetic materials. 

As we shall see later in this chapter, Eq 5.45 given in integral form below (as can be found by 
applying Gauss’ theorem to the differential form), 


B • dS = 0, 


(5.46) 


is important in determining the boundary condition on B between two materials. 


5.5.2 Mutual inductance 

A current flowing around a loop of wire (e.g. Coil 1 in Fig. 5.8) produces a magnetic field, and 
some of the field lines may pass through a second loop of wire (e.g. Coil 2). If we double the 
current I\ the magnetic field will double, as will the number of field lines through Coil 2. The 
way to quantify this threading of Coil l’s field through Coil 2 is by the mutual inductance M 
defined as the magnetic flux through Coil 2 due to current A, divided by I\ 
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M = 2^ 
h 


U 2 a i(r2) '* 2 = H i^iSr* 2 


(5.47) 


where R 21 is the separation vector from line element dri to line element dr 2 (Fig. 5.8b). Hence 


M 


Mo 

47r 


dri • dx2 
#21 


(5.48) 


and we see that the mutual inductance could be calculated also from magnetic flux due to Coil 
2 through Coil 1 divided by I 2 with exactly the same value, i.e. M — $^ 2/^1 = ^£, 1 /^ 2 - Also 
note that M only depends on the geometry (separation, shapes and sizes) of the two coils. The 
unit of inductance is the henry (H). We shall return to the topic of inductance in “Essential 
Electrodynamics” in the context of Faraday’s law. 



Coil 1 


Coil 2 



Figure 5.8: (a)-Left: randomly selected magnetic field lines due to current Ii in Coil 1, some 
of which thread Coil 2. (b)—Right: geometry for calculation of the mutual inductance of the 
two coils. 


5.6 Multipole expansion of vector potential 


For the localised current distribution J (r) the vector potential and its components are 


A(r) 


Mo 

47r 


/ 


J(Q 

R 


d 3 r', 


M r ) 


Mo 

47T 


/ 


M *0 

R 


d 3 r ', 


(5.49) 
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where i — 1,2,3 correspond to x,y,z. This is similar to Coulomb’s law for a localised charge 
distribution which we expressed as a multipole expansion (Eq. 3.5). Thus 

A l r ) = itrj Ji(r,)dV + / Ji(r/)r ' r W + '' ’ ’ (5 ' 50) 

A(r) = £^ / J(r,)dV + 4^3 y'j(r , )rT / d 3 r / + ... (5.51) 

(monopole term) (dipole term) (5.52) 


In magnetostatics a localised current distribution may be considered as being made of n closed 
circuits, T^, carrying currents, and so 


a(i)= L + 4 ^x>£*' r 

2=1 2=1 

(5.53) 

(monopole term) (dipole term). 

(5.54) 

Note, the monopole term vanishes because dr' is an exact differential, and 
dipole and higher terms. In Appendix H we show that the dipole term can 

we are left with the 

be written 

A (r)=r; r , 

(5.55) 

where the magnetic dipole moment is 


1 n r 

m = - h f r ' x dr ' ( n current loops), 

2 2=1 

(5.56) 

m= - f r / xJ(r / ) d 3 r' (volume current distribution). 

(5.57) 


Consider the current loop in the x-y plane shown in Fig. 5.9(a). The cross product of the two 
vectors r' and dr' generates a parallelogram of area |r' x rfr 7 |, with half its area corresponding 
to the area of the shaded triangle in Fig. 5.9(a) labelled dS , i.e. 


-r' x dr' — dS z. 


(5.58) 
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Hence, integrating Eq. 5.56 around this loop would give 


m = IS z 


(5.59) 


where S is the area of the loop. 



Figure 5.9: (a) Current loop in xy plane, (b) Earth’s magnetic field outside the Earth approx¬ 
imated by that of magnetic dipole m at its centre, currently misaligned by 168.5° with respect 
to the Earth’s spin angular momentum L. (c) Magnetic field of a bar magnet. 
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5.6.1 Dipole magnetic field 


The magnetic field of a magnetic dipole resembles the electric field of an electric dipole (Fig. 3.3). 
The Earth’s magnetic field outside the Earth’s surface can be approximated by that of a dipole 
of moment |m| = 8.2 x 10 22 A m 2 at the centre of the Earth (Fig. 5.9b). The dipole moment 
m is directed at an angle of about 168.5° to Earth’s spin direction, so that the “real” magnetic 
North Pole is actually 11.5° from the geographic South Pole! The magnetic field of a bar magnet 
is sketched in Fig. 5.9(c) showing the magnet’s north and south poles - by convention magnetic 
field lines leave the magnet at its north pole and enter the magnet at its south pole. 


We derive the magnetic field of a dipole here using vector identities. A derivation using index 
notation is carried out in Exercise 5-7. Taking the curl of the vector potential in Eq. 5.55, 


B(r) = x (m x 
w 4 tt V r 3 


(5.60) 


Using the product rule for the curl of a cross product, 


V x (a x b) = (V • b)a - (V • a)b + (b • V)a - (a • V)b, 


(5.61) 


we have 


v x ( m x ^) = ( v • ^) m - ( v • m )^ + (^ • v ) m - ( m • v )^> 


= 0 — 0 + 0 — (m • V)-,, 


(5.62) 

(5.63) 


where we have used V • (r/r 3 ) = 0 (r > 0) for the 1st term, and that m is a constant vector in 
the 2nd and 3rd terms on the right. Then 


B(r) = — 7^-(m • V)U 
w 47 t v 'r 3 


Mo 

47T 


r( m -V)^ + ^(m-V)r 


(5.64) 

(5.65) 
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Now, firstly 


t _v 1 ( d d d \ , Z 2 

r ( m v )^ =r [ m ^ + m y-^y +m ^ ) ( x +y 


(5.66) 


= r 


rrin 


dx y dy dz 
(—3/2)(x 2 + y 2 + z 2 ) _5//2 2x + m y (—3/2)(x 2 + y 2 + z 2 )~ b ^ 2 2y 


= + m z (-3/2)(x 2 + y 2 + z 2 ) 5/2 2z 


r (m • V)-l = — 3(m • r)r^. 


(5.67) 

(5.68) 


Secondly, 


1 . . If d d d\ ^ ^ 

-3 (m • v )r = -3 (m* — + m y — + m z —j(xx + yy + zz), 

1 , 1 
.-. -3 (m • V)r = -^m. 


(5.69) 

(5.70) 


Hence, for r > 0 the magnetic field of a magnetic dipole of moment m is 


B < r > = s 


3r(m • r) — r 2 m 


(5.71) 


This has identical form to the electric field due to electric dipole of moment p (Eq. 3.19). 


5.6.2 Torques on magnetic dipoles 


The net force on a current loop (or dipole) in a magnetic field is 


Fmag — / (j) (dr x B), 


(5.72) 


which for a uniform magnetic field is 


Fmag = I ( (p dr ) x B = 0, because (p dr = 0 


(5.73) 


Although there is no net force, a dipole does experience a torque tending to align its dipole 
moment with the magnetic field direction. To find the value of the torque consider a rectangular 
current loop, with sides of length a and c, carrying current / in a uniform magnetic field as shown 
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in the two orthogonal views in Fig. 5.10(a) and (b). Using Ampere’s force law dF = I dr x B 
we will find the force on each of the four sections. 




Figure 5.10: Geometry for calculating the torque on a rectangular current loop in the presence 
of a uniform magnetic field (a) looking in the y direction, and (b) looking in the —x direction. 

From Fig. 5.10(a), we see the top and bottom sections of the loop are perpendicular to B, and so 
Ft = —IcB x and Fb = IcB x = —Ft and these two forces cancel each other. Similarly, looking 
at Fig. 5.10(b), we see that the forces on the left and right sections are Fl = — Iasm(n — 6)B y 
and F r = la sin(7r — 0)By = — Fl, and so also cancel each other such that the net force is zero 
in agreement with Eq. 5.73. 
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From Fig. 5.10(a), we see that there will be a torque on the current loop due to forces Ft and 
F# equal to N = — |Fx|asin0y = — IacB sin 0y, and this tends to align m parallel to B. Given 
that B = Bz and m — Iac(cos0z + sin0x) we find that 

N = m x B. (5.74) 


Such current loops in a magnetic field form the basic design of the galvanometer at the heart of 
all analogue current, voltage or charge measuring devices (Fig. 5.11), and also electric motors. 
Similarly, these torques also tend to align electron magnetic moments (associated with spin) 
parallel to B causing paramagnetism , to be discussed in the next chapter. 



Figure 5.11: A simple galvanometer to measure current / has coils through which the current 
flows, a permanent magnet with shaped pieces to provide the same field configuration with 
respect to the coils at all deflection angles, a spring to provide a torque to balance the torque 
on the magnetic moment of the current through the coils, and a needle pointer. 

To rotate the dipole moment, work must be done against the torque, and so the potential energy 
of a dipole in a magnetic field is 


r0 rO 

W{0) = / N(Q')d6' — mB / sin d'dO' = m£(cos0 ref - costf), (5.75) 

J 0-rp f " $ref 

(5.76) 


' ^ref 

W = - m B 


where, for convenience, we have chosen the reference angle to be d Te f = n/2. 


As an aside, although for a uniform magnetic field there is no net force on a magnetic dipole, 
in a non-uniform magnetic field there is. The force on the dipole is 


F = 


—VVF = V(m-B). 


(5.77) 
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If the dipole is aligned with B, it will experience a force in the direction of increasing B. 
Actually, the result in Eq. 5.77 glosses over a number of issues including whether the dipole 
is due to a current loop, electron orbital motion or spin. Nevertheless, Eq. 5.77 explains why 
paramagnets and ferromagnets are attracted to strong magnetic fields whereas diamagnets are 
repelled. 

Summary of Chapter 5 

Magnetic forces 

— Lorentz force law F = Q(E + vxB). 

— Force on a circuit in a magnetic field F mag = ^(/drxB). 

— For surface and volume currents F mag = J K x B = / J x B d 3 r. 

Magnetic field and vector potential 

— Biot-Savart law 



(5.78) 


Vector potential A where B = V x A 



(5.79) 


Ampere’s law and the no magnetic charge law 


Ampere’s law 




Magnetic flux = J B(r) • dS = j) A(r) • dr 


No magnetic charge law V • B = 0 
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Multipole expansion of vector potential 


— The monopole term is zero, and the dipole term is 


A(r) = 


/io m x r 
47t r 3 


(5.81) 


Magnetic dipole moment for n current loops and for a volume current distribution 


m 


= ^ I{ (j) r ' x dr m = ^ f r ' x J(r ') cfir'. 

i=l '' r * ^ 


(5.82) 


— Dipole magnetic field 


B « - s 


3r(m • r) — r 2 m 


(5.83) 


Torque on dipole N = m x B. 
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Exercises on Chapter 5 

5-1 Two parallel wires are separated by distance a and carry currents I\ and I 2 in the same 
direction. Find the force per unit length of wire. Include a diagram showing the direction 
of the force. If I\ = I 2 = 1 A and a — 1 m, what is the magnitude of the force per unit 
length? 

5-2 Using the equation for the vector potential in terms of the current density (Eq. 5.21) find 
the vector potential of an infinite straight wire along the 2 axis carrying current I. 

5-3 We can add the gradient of a scalar field U( r) to A(r) without changing B(r). This is 
called a gauge transformation. The “gauge” of the vector potential is determined by the 
value of V • A. Show that in magnetostatics V • A = V 2 U(r). 

5-4 (a) Find the vector potential of the constant magnetic field Bq = (b) check 

that the vector potential you find does give the desired magnetic field, (c) find V • A and 
check it is what is expected in magnetostatics. [Hint: first find the vector potential A^(r) 
of the simpler constant field B® z by writing down the components of V x A^ in Cartesian 
coordinates before appealing to the symmetry of the problem, and then integrating.] 

5-5 A steady current / flows down a long cylindrical wire of radius b. Find the magnetic field 
both inside and outside the wire. 

5-6 A semi-infinite solenoid of radius a, has n turns per unit length, extends from z = —00 to 
z = 0 along the 2 : axis and carries current / in the +<fi direction. Magnetic flux is confined 
to the solenoid, but emerges isotropically from its end at the origin as shown below. 



(a) On the cone of half-angle 6 with apex at the origin there is a circular loop (as shown) 
with all points on the loop being at distance r a from the origin. Find the magnetic flux 
passing through this loop, (b) Find the magnetic vector potential at the point (r, #,</>), 
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and take it’s curl to find the magnetic field. [The expression for magnetic flux through a 
loop in terms of the vector potential may be useful here.] 

5-7 By taking the curl of the vector potential for a magnetic dipole with moment m located 
at the origin, find it’s magnetic field using index notation. 

5-8 A circular current loop in the xy plane has radius a and is centred on the origin. It carries 
current / in the ^-direction. There is a uniform magnetic field B(r) = £>o(cos#z + sindy) 
present. By integrating the torque dN = r' x dF on line element dr' of the current loop 
at r', find the torque on the entire current loop. Compare your result with the result you 
would get by first finding the current loop’s dipole moment, and then applying the formula 
for the torque on a magnetic dipole in a magnetic field. 
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6 Magnetism of materials 

Learning objectives 

— To learn that a material placed in a magnetic field acquires a magnetic dipole moment 
per unit volume described by the magnetisation field, and that magnetisation results in 
surface and volume magnetisation currents. 

— To be able to calculate the magnetisation currents, magnetic field and the magnetic dipole 
moment of a permanently magnetised object. 

— To understand that the magnetisation currents must be added to any free currents present 
when applying Ampere’s law, and that this leads to the introduction of a new field H. 

— To be able to apply the boundary conditions on B and H at an interface between two 
materials. 

— To learn that for many materials the magnetisation is proportional to H, with the con¬ 
stant of proportionality being the magnetic susceptibility, which in turn is related to the 
permeability. 

— To understand the causes of paramagnetism, and how a semi-classical treatment leads to 
Curie’s law. 

— To understand the causes of ferromagnetism, and that it requires no applied magnetic 
field to sustain it. 

— To know how the magnetisation of a ferromagnetic material changes with applied magnetic 
field or H, and to be able to calculate the work done as a ferromagnetic material is brought 
around a hysteresis loop. 

— To understand the concept of magnetic circuits, and to be able to calculate magnetic fields 
in simple magnetic circuits. 


6.1 Magnetic materials 


While in an external magnetic field, a paramagnetic substance becomes magnetised with a 
net dipole moment parallel to the magnetic field. Paramagnetism occurs if, for a particular 
substance, the atomic electrons have unpaired spin allowing the electron magnetic moments 
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of those unpaired electrons to experience a torque aligning them with B. Ferromagnetism is 
due to the quantum mechanical exchange interaction causing long-range ordering of electron 
spins within a “domain” — a concentration of aligned magnetic dipoles. If a material is neither 
paramagnetic nor ferromagnetic, it will show (usually weak) diamagnetism. This is due to the 
effect of the magnetic field acting on orbiting electrons and changing the magnetic moment 
associated with their orbital motion in such a way that the net dipole moment becomes anti¬ 
parallel with B. Diamagnetism occurs for all atoms, but in paramagnetic and ferromagnetic 
materials is masked by the typically stronger paramagnetism and ferromagnetism. 


6.2 Magnetisation field 


In the presence of a magnetic field matter becomes magnetised in an analogous way to that in 
which matter becomes polarised in the presence of an electric field. Magnetic dipole moments 
due to electron spin or orbital motion get partially or fully lined up in some direction giving 
the material a net magnetic dipole moment. To describe this we define the magnetisation field 
M(r) in a similar way to the polarisation field in electrostatics, namely the magnetisation field 
is defined as the magnetic dipole moment per unit volume. The magnetisation field M has 
units of magnetic dipole moment (A m 2 ) divided by those of volume (m 3 ), i.e. its units are (A 
m -1 ). In the same way that in electrostatics the electric field of a polarised object is identical 
to that due to its polarisation charges, the magnetic field of a magnetised object is identical to 
that of its magnetisation currents , also known as “bound currents”. 


6.3 Magnetisation currents and Ampere's law 


The vector potentials of a magnetic dipole with moment m at the origin, of a magnetic dipole 
with moment m located at r', and of a magnetised object are given by 


A(r) 


fi o m x r 
47r r 3 ’ 


A(r) 


fi o m x R 
An R 3 ’ 


A(r) 


/fo f M(rQ x R 3 , 

47 T J i? 3 


(6.1) 


where, as usual, R = r — r'. 


To derive the magnetisation currents of a magnetised object we will need the following: an 
identity for V'i? -1 , the product rule for V x (aF), and a corollary to Gauss’ Theorem, all 
of which are obtained in Exercise 6—1. The derivation of magnetisation currents proceeds as 
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follows. First re-write the cross product in the integral in Eq. 6.1 


A ( r ) = £/ M ( r/ ) x 


( 6 . 2 ) 


Then, use V'i? 1 = +R 2 R (Eq. B.9) to find 


A(r)= S/ M(r ' ) 


xv'fl WV. 


(6.3) 


Next, using the rule for the curl of the product of scalar and vector fields (Eq. 1.9) we obtain 


A(r) = ^ 

47 T Jy R 


L- 


xM ^/ /* f 


(f r - ^ 

47T J V 


/ v ' 


R ) 


(6.4) 


Finally, using the corollary to Gauss’ theorem f v V x Fd 3 r = — j> s F x dS we get 


A(r) = 


Mo 

47T 


r v: 

JV 


A(r) ^£X- 


+ ^0 / xdg/ 

^ t Js\R ) 

x M d 3 / . /fo / MxdS ; 
i? 47T /X -R 


(6.5) 

(6.6) 


Hence, the field of a magnetised object is the same as would be produced by “magnetisation 
currents” 


Jmag(r) = V x M(r) (volume), K mag (r) = M(r) x n (surface). 


(6.7) 


To visualise these magnetisation currents, consider a rectangular block of material uniformly 
magnetised in the —x direction as in Fig. 6.1(a). We can imagine the magnetisation field as 
being due to the tiny current loops shown. If M is constant throughout the material then each 
of these fictitious current loops carries an identical current, such that inside the material the 
current on one side of a loop is in the opposite direction to the current on the nearest side of an 
adjacent loop and cancels it out (see inset in Fig. 6.1a). So, J m ag(r) = 0 inside the uniformly 
magnetised material, as it must be since V x M = 0. 

There is no magnetisation surface current density in the surfaces perpendicular to M (with 
darker shading), i.e. parallel to the yz plane in Fig. 6.1(a), because M x n = 0 there. However, 
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for the loops which have one side next to the surface, the currents on that side of the loops have 
no opposing currents to cancel with. The result is a magnetisation surface current K mag (r) = 
M(r) x n such that on the surface of the side with lighter shading parallel to the xz plane 
K mag (r) — M z whereas on the top surface with lighter shading (parallel to the xy plane) 
K m ag(r) = My. 



Figure 6.1: Physical interpretation of magnetisation currents, (a) Uniform magnetisation. 
Inset shows detail of cancelation of internal currents (b) M x varies with y showing: left - 
internal currents don’t cancel; right - the relative magnitudes of the surface magnetisation 
current densities K mag as it varies over the surfaces, as well as the direction of the volume 
magnetisation current density J mag . 
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If M varies with position inside the material, as well as the surface magnetisation currents 
there will also be a volume magnetisation current density provided V x M(r) ^ 0. One 
example is is shown in Fig. 6.1(b) where M(r) depends only on y and varies linearly with ?/, 
i.e. M(r) M x (y) x — (Mq + ay) x such that J(r) = V x M(r) —dM x /dy z a z . The 
resulting surface and volume current densities are illustrated in Fig. 6.1(c) where on the top 
surface K mag (r) = M x (y) y which varies with position, while on the sides parallel to the x-z 
plane the magnetisation surface current densities are constant with the values K mag (x,0, z) — 
—M x ( 0) z Mo z and K mag (x, d, z) — M x (d) z —(Mg + ad) z. 


6.3.1 Cylindrical magnets 


As an example of magnetisation currents, we consider a permanent magnet in the form of a 
short cylinder with uniform magnetisation M along its axis as in Fig. 6.2(a). The magnetisation 
surface current K mag (r) = M(r) x n is zero on the two ends, but on the cylindrical surface 
Kmag(r) = M<p (Fig. 6.2b). This surface magnetisation current distribution would give an 
identical magnetic field to a free surface current density K — nl (p where I — M/n is the 
current in a tightly-wound solenoid of the same size as the magnet (Fig. 6.2b), and n is the 
number of turns per metre. 

(a) (b) 




Figure 6.2: Magnetic field of (a) a short cylinder with uniform magnetisation M = Mz parallel 
to its axis, and (b) a solenoid of identical size with n coils per unit length and carrying current 
I = M/n, effectively a surface current K = nlcfi. 


116 


Download free eBooks at bookboon.com 














































Essential Electromagnetism 


Magnetism of materials 


6.4 Ampere's law in magnetised materials 


For magnetic materials we must include J mag as well as the “free current” Jj, i.e. J used 
previously must be replaced by the total current density as free currents and magnetisation 
currents are equally good at producing magnetic fields. Hence, 


V x B = /i 0 (J / + Jmag ) — fJLoQ f + V X M). (6-8) 

Hence, Ampere’s law in magnetised materials is 

VxH = J/, where H = — - M. (6.9) 

tM) 

H is sometimes called magnetic intensity but more often simply “H”. It has units Am 1 . Using 
Stokes’ theorem we obtain Ampere’s law in integral form 

j) H • fir = If. enc i, (6.10) 

where If , enc i is the free current enclosed by loop T. 


6.4.1 Boundary Conditions on B and H 


We obtain these boundary conditions from the no-magnetic-charge law and Ampere’s law in 
integral form. Consider the Gaussian pillbox shown in Fig. 6.3(a) but with infinitesimally 
narrow sides. The surface integral of B over the pillbox with identical top and bottom surface 
areas S± = S 2 is given by 


j) B(r) • dS = 0, (Bi • niSi + B 2 • r^S^) = 0. (6.11) 

B above = B below> ( 6 - 12 ) 

and we see that the component of B normal to the interface is unchanged across the boundary. 

To find the boundary condition on H we apply Ampere’s law f H • dr = If,end 1° Amperian 
loops r a and of length d and infinitesimal width (Fig. 6.3b). Consider loop T a whose sides 
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Figure 6.3: Geometry for obtaining the boundary conditions on (a) the normal component 
of B, and (b) the component of H parallel to the interface between material 1 (above the 
interface) and material 2 (below). The thickness of the Gaussian pill box and the widths of the 
Amperian loops are exaggerated in the diagrams, and should be considered as infinitesimal. 


are parallel to Kj such that enclosed current is zero, and loop r& whose sides are perpendicular 
to K f and has enclosed current Kf d, 


(f) H • dr = (Hi -t a d — H 2 • t a d) = 0, 
Jr a 


(6.13) 
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and 


(f H • dr = (Hi -t b d — H 2 • t b d) = K f d. (6.14) 

Jr b 

So the component of H parallel to both the surface and Kj is unchanged, and the component 
of H which is both parallel to the surface and perpendicular to K f is changed by an amount 
Kf. Combining the results for the two parallel components of H we obtain 

( H Lve- H Lw ) = K/xn. (6.15) 


6.5 Magnetic susceptibility and permeability 

For linear materials, M = \M H where \M is the magnetic volume susceptibility , and so 


B — /i-o(H + M) — /+)(! + Xm)H — /iH 


(6.16) 


where fi is the permeability of the material, and fi o is the permeability of free space. Since the 
magnetic dipoles are at the atomic/molecular level, with M depending on the number density 
of atoms/molecules (which varies by orders of magnitude between gases and solids), it makes 
sense when comparing the magnetic properties of materials to divide the magnetic volume 
susceptibility by the density. This quantity, called the magnetic mass susceptibility is tabulated 
in Table 6.1 for several elements, compounds and substances. 


Substance 

Xm/p 

Substance 

Xm/p 

Substance 

Xm/p 

Substance 

Xm/p 

Aluminium 

+0.82 

Bismuth 

-1.70 

Copper 

-0.107 

Gallium 

-0.30 

Germanium 

-0.15 

Gold 

-0.19 

Helium 

-0.59 

Hydrogen 

-2.49 

Krypton 

-0.41 

Lead 

-0.15 

Mercury 

-0.21 

Neon 

-0.41 

Nitrogen 

-0.54 

Oxygen 

+133.6 

Phosphorus 

-1.13 

Platinum 

+1.22 

Potassium 

+0.65 

Silicon 

-0.16 

Silver 

-0.25 

Sodium 

+0.75 

Sulphur 

-0.62 

Uranium 

+2.19 

C0 2 

-0.59 

H 2 0 

-0.90 

HC1 

-0.75 

H 2 S0 4 

-0.50 

NaCl 

-0.64 

nh 3 

-1.38 

NO 

+59.3 

Perspex 

-0.5 

Polyethylene 

+0.2 

P.V.C. 

-0.75 


Table 6.1: Mass susceptibilities, Xm/p (10 -8 kg -1 m 3 ) at 20° C of some common elements, 
compounds and materials (data taken from Kaye & Laby 
http://www.kayelaby.npl.co.uk/general_physics/2_6/2_6_6.html). 
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6.6 Paramagnetism 


Pierre Curie discovered that for ordinary temperatures and magnetic fields paramagnetic ma¬ 
terials obey xm = (M/H) oc 1/T. Paramagnetism is analogous to orientational polarisability 
of permanent molecular electric dipoles. Magnetic torques on atomic dipole moments cause 
paramagnetism, and it is the magnetic moments associated with unpaired spin that are respon¬ 
sible. Unpaired electron spin will occur where there is an unfilled atomic sub-shell. However, 
in metals only those electrons with energies near the Fermi energy can respond in this way and 
so paramagnetism in metals is weak and does not follow the Curie law. 


For a particle of mass M p and charge q in a circular orbit of radius R with speed v (Fig. 5.4) 
the magnitudes of the magnetic dipole moment mdip and angular momentum L are 


?7idip — and L — M p vR , 


(6.17) 


where the dipole moment is the product of the current (q divided by orbital period 2ttR/v) and 
the area of the orbit ( nR 2 ). The gyromagnetic ratio is defined as 


^dip _ Q 


2M p 


and 777/dip = 


qL_ 

2 Mr, 


(6.18) 


For atoms we need to take account of quantum mechanics, in which the quantum of angular 
momentum is h=h/2n where h is Planck’s constant. The natural unit of 777/dip is therefore the 
“Bohr magneton” 


M B 


e h 
2m e 


(6.19) 


The total angular momentum J = L + S is the sum of the orbital angular momentum of electrons 
and their spin. For B = Bz the quantum number M for the ^-component of J is relevant here, 
such that the z-component of the magnetic dipole moment mdi p is 


^dip,2; 


-gj/LL B M, 


gj = i + 


J {J + 1) + S(S + 1) — L(Ij + 1) 
2J(J+1) 


( 6 . 20 ) 


where gj is the “Lande g- fact or”. Note that if L = 0 (no orbital angular momentum) then 
J = S and gj — 2 for electron spin, and that if S = 0 (no spin angular momentum) then J — L 
and gj — 1 for electron orbital angular momentum. 
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As an aside, the Zeeman effect is the atomic energy-level splitting in a magnetic field 

A E = -m dip • B = +gjn B MB (6.21) 

and this provides an important method for radio astronomers to determine the magnetic field 
present in some astrophysical sources. 


6.6.1 Derivation of susceptibility 

Paramagnetism is due to partial alignment of permanent atomic magnetic dipoles. For a full 
treatment, a quantum statistical mechanical approach is required, taking account of the quan¬ 
tum numbers for total angular momentum J, orbital angular momentum L, and spin S. How¬ 
ever, the problem is simplified for the case atoms with zero orbital angular momentum, as we 
then have gj — 2 and only two possible spin states M s = ±^. The corresponding magnetic 
dipole moments are 



MAERSK 


I joined MITAS because 
I wanted real responsibility 


The Graduate Programme 
for Engineers and Geoscientists 

www.discovermitas.com 


I was a construction 


supervisor m 
the North Sea 
advising and 
helping foremen 
solve problems 


Real work 
International opportunities 
Three work placements 


Sir 


50 


0 


121 

Download free eBooks at bookboon.com 











Essential Electromagnetism 


Magnetism of materials 


^di P ,z = ~9j^bM s = -2 fi B ( ±- ) = T 9b, 


( 6 . 22 ) 


and the corresponding potential energy values are 


W = -m • B = -m d[ ^ z B = ±ii b B. 


(6.23) 


Using the Boltzmann distribution, for temperature T the probabilities of the electron magnetic 
moment being parallel and antiparallel to the magnetic field are 


^parallel = A e^ B / kT , Parallel = A e~^ B / kT 


(6.24) 


where Po is a constant and (Pparaiiei + ^antiparaiiei) — 1- Then, for number density N of unpaired 
electrons the magnetisation is 


(-^parallel T^ntiparallel) N flB B 


e HBB/kT _ e —fi B B/kT 
gUsB/kT g—UBB/kT 


N fi b B. 


(6.25) 


For i±bB <C fcT, i.e. for ordinary temperatures and magnetic fields, 


M ~ H), 


(6.26) 


and the magnetic susceptibility becomes 


Xmn 


NfiQli 2 B 

kT ’ 


(6.27) 


consistent with the Curie law. 


6.7 Ferromagnetism 


Ferromagnetism is non-linear , and no external fields are needed to keep a ferromagnet magne¬ 
tised. Ferromagnetism is due to long-range ordering of unpaired electron spins causing then to 
line up parallel to each other. Anti-ferromagnet ism is similar, except that neighbouring dipoles 
point in opposite directions. Ferrimagnetism has oppositely pointing neighbouring dipoles of 
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different magnitude which don’t cancel. These three possibilities are illustrated schematically 
in Fig. 6.4. 


Ferromagnetic Antiferromagnetic Ferrimagnetic 


cjp cjp cjp cjp cjp cjp 

cj? <^p cjp <^p c|p ijp 

cjp ^ cjp ^ cj 

b 

cjp <jp cjp cjp cjp cjp 

cjp cj? <^P <^p 

cjp ^ cjp ^ cj 

b 

cjp cjp cjp cjp cjp cjp 

<j> cjp c|p cjp cjp cjp 

4*+4 + *\ 

b 


Figure 6.4: Arrangement of dipoles in crystal lattice planes of ferromagnetic, antiferromagnetic 
and ferrimagnetic materials. 

Alignment occurs in microscopic “domains” with the magnetisation direction of each domain 
usually being random, and giving an overall unmagnetised state. When placed in an external 
magnetic field, domains with magnetisation aligned partially with B grow at the expense of less 
favourably aligned ones so that the average magnetisation is in the direction of B. Torques N = 
m x B try to align individual dipoles with B, but they like to stay parallel to their neighbours, 
and a lot of work is needed to shift the domain boundaries to achieve greater alignment. For 
ferromagnetic materials, the magnetic field B m = H oM due to the resulting magnetisation M is 
typically much larger than the applied field Bq = /ioH. In ferromagnetism and ferrimagnetism 
the material becomes paramagnetic above a critical temperature at which a phase transition 
occurs called the Curie temperature (after Pierre Curie) and in anti-ferromagnetism above the 
Neel temperature (after French physicist Louis Neel 1904-2000). Practical applications of this 
phase transition occur in electro-optical devices and hard disk drives when heat is used to erase 
and write data. 

Typically one would examine the magnetic properties of a ferromagnetic material by placing 
the sample inside a solenoid and varying the current (and consequently H) sinusoidally, often 
at mains frequency (50 Hz or 60 Hz), and displaying Bm vs. H on an oscilloscope or computer 
screen. An example is shown in Fig. 6.5(a) for a magnetically-soft iron-based amorphous alloy. 
What is found is that Bm remains positive even when H returns to zero following its maximum 
positive value, and that Bm reaches zero only after a time lag when H is sufficiently negative. 
Provided iL max is not too large Bm vs. H is approximately linear. Such a material would be 
used mainly for magnetic cores in transformers. For larger applied fields, M saturates when all 
domains having spins partially aligned with H have grown as much as is possible. 

This type of non-reversible process is called hysteresis. To bring the sample through one cycle 
of the “hysteresis loop” work must be done per unit volume of the ferromagnetic material equal 
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H (kA m' 1 ) H (kA m' 1 ) 


Figure 6.5: Hysteresis in (a) a magnetically-soft iron-based amorphous alloy, (b) a magnetically- 
hard alloy of iron, aluminium, nickel and cobalt, Alnico. 


to the loop’s area, W = j) H dB. Depending on the practical application of the magnetic 
material, which can vary enormously from computer discs to permanent magnets to magnetic 
films for fridge magnets to magnetic shielding of electronics, one would choose a material which 
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has “low hysteresis” (area of hysteresis loop is small), and is termed “magnetically soft” (good 
for magnetic cores for transformers and for electromagnets), or one which has “high hysteresis” 
and is termed “magnetically hard” (good for permanent magnets), or one which has the right 
shape of the hysteresis loop for the job in question. 

The hysteresis loop of a magnetic material is characterised by its remanence Br (T), its coerciv- 
ity He (A m -1 ) and its saturation Bs (T), and these three quantities are labelled in Fig. 6.5(a). 
For this particular amorphous alloy: Br = 0.1 T, He — 10 A m -1 and B$ — 1 T. Again, de¬ 
pending on the application, the shape of the hysteresis loop is important - it depends not only 
on the material but also on whether it has been annealed, i.e. heated and cooled in the presence 
of a magnetic field aligned in a particular direction. An example of a hysteresis loop with a 
very different shape to that of the magnetically-soft amorphous alloy is shown in Fig. 6.5(b), 
in this case for a magnetically hard alloy of iron, aluminium, nickel and cobalt ”Alnico”, which 
would be used for permanent magnets. 

A recent development allows us to image ferromagnetic domains using scanning electron micro¬ 
scopes (SEM). These microscopes work by having a beam of keV to multi-keV energy electrons 
scan across the surface of the electrically conducting sample in a raster scan pattern. If the 
object is non-conducting it must be coated with a thin layer of a good conductor such as gold. 
Some of the scanning electrons are absorbed and the resulting current depends on the local 
surface height and material properties, and is used to construct an image of the surface which 
is displayed on a computer screen. Electrons scattered off the surface can also be collected in 
real-time, as can X-rays emitted as a result of the surface bombardment, and are also used in 
construction of the image in some SEMs. 

Recently, a detector was developed to collect and analyse the spin state of the scattered elec¬ 
trons, which in turn depends on the magnetisation direction of the patch of surface being 
imaged. This scanning electron microscope with polarisation analysis (SEMPA) technique en¬ 
ables images of the domain structure on the surface of a ferromagnetic material to be made. 
Figure 6.6 shows images of magnetic domains and domain wall dynamics in Metglas® ribbons. 
SEMPA measurements of the magnetic domains are shown at four points on the hysteresis loop 
(centre) showing the magnetic state corresponding to each image. 
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Figure 6.6: Hysteresis in an amorphous alloy. Red (darker shading) indicates posi¬ 
tive magnetisation, green (lighter shading) indicates negative. The small, relatively im¬ 
mobile domains are pinned by stress fields induced during fabrication. Each of the im¬ 
ages is approximately 980 microns wide. Image is courtesy of A. Gavrin et al., from 
http://www.nist.gov/cnst/epg/wall_motion_proj.cfm (reproduced with permission). 


6.7.1 Magnetic circuits 


Consider two rods made of plastic, each of which is heated in order to bend it into a horseshoe 
shape. The two plastic horseshoes are placed end-to-end almost touching but leaving two short 
air gaps between them. Wires are wrapped around part of one of the horseshoes to form a 
short solenoid as shown in Fig. 6.7(a). A current passing through the solenoid will produce a 
magnetic field which, because the magnetic susceptibility of plastic is tiny, will be only slightly 
different to the field the solenoid would produce in air. 



Figure 6.7: (a) Two horseshoes made of plastic and placed end-to-end almost touching but 
leaving small air gaps, and having a short solenoid wound around part of one horseshoe, (b) As 
part (a) but with the horseshoes being made of ferrite core material. 


Now consider horseshoes of identical size and shape to the plastic horseshoes, but made of 
ferrite core material. Ferrite is magnetically soft and has fi /jlq. There is an identical short 
solenoid with N coils carrying current / as shown in Fig. 6.7(b). The magnetic field in the 
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ferrite core will be dramatically different to that in the plastic horseshoes as can be shown 
using the boundary condition on H at the boundary between the ferrite core material and the 
air. Since there are no free surface currents in the ferrite, the component of H parallel to the 
surface, HJI, must be continuous across the boundary between ferrite and air 

® ferrite/ /ferrite = ®Iir/Mair- (6.28) 

Hence, b| rrite > Bj ir as /iferrite ^ l 1 o> and the magnetic field will be largely confined to the 
ferrite core, with the field lines making loops as shown. Using the boundary condition on the 
normal component of B, the magnetic field in the two air gaps between the two horseshoes 
must be unchanged across the boundaries as the air gap surfaces are normal to B inside the 
ferrite core. 

As the ferrite core effectively channels the magnetic field, effectively we have a magnetic circuit . 
As an example, we shall calculate the magnetic field in the air gaps for a permeability of 
ji i = 5,OOO^o for the horseshoe on the left and (12 = 3,000/io for the horseshoe on the right - 
strictly speaking ferromagnetic materials are non-linear, but we shall assume these values apply 
to initial magnetisation before the field saturates. The other values we shall use in this example 
are / = 0.2 A, N — 50, the lengths around the two horseshoes being L\ — 30 cm and L 2 = 20 
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cm, a < L 2 and the air gaps having width d — 0.2 mm (<C a). Applying Ampere’s law to the 
closed curve following the centre of the cross-sections of the ferrite core horseshoes 


B B B 

H • dv — /gnci 1 .*. dj\ — Z/2— 2 d — — TV I . 

Ml M2 Mo 


Substituting the numbers we find B = 0.024 T. 


(6.29) 


Summary of Chapter 6 


Macroscopic theory of magnetism 

— The magnetisation field M is defined as the magnetic dipole moment per unit volume. 

— Paramagnets have M parallel to the applied field, diamagnets have M anti-parallel. 

— Magnetised materials have magnetisation currents 

Jmag(r) = V x M(r) (volume), K mag (r) = M(r) x n (surface). (6.30) 

— Ampere’s law in magnetised materials V x H = Jj, ^ H • dr = //, enc i, 
where H = B/ /io — M. 

— Boundary Conditions between magnetic materials 

B above = B below> (Hj bove - Hj elow ) =K/Xn, (6.31) 

— Linear materials: susceptibility \M and permeability fi 

B — /io(H + M) = /i o(l + xm)H = /iH. (6.32) 


Paramagnetism 

— Magnetic torques on atomic dipole moments cause paramagnetism. 

— At normal temperatures and applied fields partial alignment occurs giving 


XM = 


B 


(±§nmn 
3 kT 


(6.33) 


128 


Download free eBooks at bookboon.com 




Essential Electromagnetism 


Magnetism of materials 


Ferromagnetism 

— Ferromagnets are non-linear due to long-range ordering causing unpaired electron spins 
to line up parallel with each other. 

— Alignment occurs in microscopic domains with the direction of M in each domain being 
random, and giving overall unmagnetised state. 

— When a magnetic field is applied, torques N = mxB try to align individual dipoles with 
the applied field. 

— Work is done to shift the domain boundaries to cause overall magnetisation. 

— This non-reversible process is called hysteresis - a plot of Bm vs. applied H shows the 
hysteresis loop. 

— The area of the loop gives the work done to take the sample around one cycle of the loop. 


Exercises on Chapter 6 


6-1 Derive the following which are needed to obtain the magnetisation currents of a magnetised 
object: (a) Identity for V'i? -1 , 

v'ir 1 = +ir 2 R. (6.34) 

(b) Product rule for V x (aF), 

V x (aF) = (Va) x F + aV x F. (6.35) 


(c) Corollary to Gauss’ Theorem, 


L 


V x F d 3 r 


|fx dS. 


(6.36) 


6-2 A thin disc of magnetised material is coincident with the xy plane. It is of thickness s and 
radius a and has magnetisation M = Mqz. Find the magnetisation current, and from this 
find the magnetic dipole moment of the disc. Compare this with what you would get by 
multiplying the disc’s volume by M. 
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6-3 Consider a permanent magnet in the form of a short cylinder of radius a extending along 
the £ axis from z — —L to z — +L and having uniform magnetization M = Mqz. (a) Find 
B and H at all points (0,0, z) on the cylinder’s axis, and plot B(0,0,z) and H(0,0,;z) vs. 
z. (b) Discuss whether the result obtained in part (a) obeys Ampere’s law for H in integral 
form. 

6-4 A cylindrical rod of radius a and length h a is permanently magnetised along its length 
which coincides with the z direction, i.e. M = Mqz. (a) Find the surface magnetisation 
current K mag (r), and use it together with Ampere’s law to find B and H inside and outside 
the rod (assume h —oo). (b) The rod is now bent into a circle of circumference (h + 2 L) 
such that there is an air gap of width 2 L < a. Plot B and H along the axis of the magnet 
in the vicinity of the air gap for the case of a — 0.5 and L = 0.2. 

6-5 Consider a permanently magnetised sphere of radius a with uniform magnetisation M(r) = 
Mqz. (a) Find the surface magnetisation current density, and use this to find the magnetic 
dipole moment of the sphere. Compare this with what you expect given the volume of the 
sphere and the magnetisation field, (b) Find B and H at the centre of the sphere. 

6-6 Consider the hysteresis loops of the magnetically-soft iron-based amorphous alloy and the 
magnetically-hard alloy of iron, aluminium, nickel and cobalt shown in Fig. 6.5. 

(a) Estimate the work done to bring 1 cm 3 of each material through one cycle of the 
hysteresis loop, (b) Two transformers operating at 50 Hz have magnetic cores of volume 
100 cm 3 (one of each type of material) and are (unwisely) operated at a current at which 
saturation occurs. How much power is lost as heat in each case? [You could print Fig. 6.5 
and estimate the area by drawing over it a grid and measuring by hand sufficient points 
on the graph to get within say 20% accuracy for the area.] 
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Due to limited space available in this series, it is not possible to go into the material in great 
depth, so I have attempted to encapsulate what I consider to be the essentials of electrostatics 
and magnetism. In a separate volume to be published by Ventus Publishing ApS, “Essential 
Electrodynamics”, I cover in similar depth Maxwell’s equations, electromagnetic waves and 
radiation. Both books contain student exercises, and I give the solutions to these exercises in 
separate volumes also to be published by Ventus: “Essential Electromagnetism - Solutions” and 
“Essential Electrodynamics - Solutions”. 

My aim with the present book is not to replace existing textbooks in these subjects of which 
there are many excellent examples, but rather to supplement the textbook specified by the 
professor teaching the course — I have always found the need to read about a particular topic in 
at least two books to improve my understanding. Nevertheless, if appropriately supplemented, 
this book together with “Essential Electrodynamics” could serve as a textbook for 2nd and 3rd 
year electromagnetism courses at Australian and British universities, or for junior/senior level 
electromagnetism courses at American universities/colleges. 

A few of the books that I found helpful over many years of teaching electromagnetism are listed 
below. 


— Barger, VD & Olsson MG 1987 Classical electricity and magnetism: a contemporary 
perspective , Allyn and Bacon, Boston 

— Bleaney, BI & Bleaney, B 1976, Electricity and magnetism , 3rd Edn. Oxford University 
Press, Oxford 

— Cheng, DK 1989, Field and wave electromagnetics , 2nd Edn. Addison-Wesley, Boston 

— Feynman, RP, Leighton, RB & Sands, M 1964, The Feynman Lectures on Physics, Volume 
77, Addison-Wesley, Reading, Mass. 

— Griffiths, DJ 1981, Introduction to Electrodynamics , 3rd Edn. Prentice Hall, New Jersey. 

— Heald, MA & Marion, JB 1994, Classical Electromagnetic Radiation , 3rd Edn. Brooks/Cole, 
Pacific Grove CA 

— Jackson, JD 1998, Classical Electrodynamics , 3rd Edn. Wiley, New York 
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A Vectors 


Vectors have both a magnitude and a direction. In print they are represented by a letter in 
bolface type, e.g. the magnetisation might be written M, or by an arrow over the symbol when 
handwritten, e.g. Ai. In three-dimensional space we can have vectors such as the velocity of a 
particle at time £, say v(t), and vector fields such as the flow velocity within a fluid v(x, ?/, z, t) 
at position (x, y, z) at time t. One can also have scalars such as the charge q of a particle, and 
scalar fields such as charge density p(x,y,z,t) at position (x,y,z) at time t. What scalars and 
scalar fields have in common is that they have just one component (it takes only one number 
to define them), whereas both vectors and vector fields have a magnitude and a direction, and 
so it takes three numbers to define them. 

Perhaps the simplest vector to describe is the position vector r of some point in space defined 
by its Cartesian coordinates (x,y,z), or the vector r' being the position vector corresponding 
to the point (a/, y' , z f ). The position vector r' is represented by a an arrow extending from the 
origin (0,0,0) to (x f ,y f ,z f ) as in Fig. A.1(a). Its magnitude is just the length (measured in m) 
of the arrow which, using Pythagoras’ theorem, is r' = \Jx' 2 + y' 2 + z' 2 . The vector’s direction 
is represented by three direction cosines , cos 6^, cos 6' y and cos 6^, which are the cosines of the 
angles the vector makes with the three Cartesian axes as shown in Fig. A. 1(b). 


(a) (b) (c) 



Figure A.l: (a) Position vector r' of point (x',y',z'). (b) Position vector r' and the angles 
6 X1 6'y and 0' z \ note that x' = r' cosO x , y' = r' cos0' y and z' = r' cos 6^. (c) Unit vectors in the 
x, y and z directions. 

We define a unit vector as being a vector having a dimensionless magnitude of 1. We have 
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just seen one example — the direction cosines of r' are the Cartesian components of the unit 
vector in the direction of r', i.e. r' = (cos 0' x , cos 0' y , cos 9' z ). Note that the “hat” above vector r' 
indicates that it is a unit vector, i.e. in this case 



OVy'V) 


V ®' 2 + V 2 + ^' 2 


(cos 6' x , cos O' , cos 9' z ). 


(A.l) 


A unit vector r' pointing in the direction of r' has been added to Fig. A. 1(b). Unit vectors in 
the x, y and z directions are written x, y and z (in some books they are written i, j and k). 
These three unit vectors form an orthogonal set of basis vectors (Fig. A.lc). 


The force F on a particle has a magnitude F (measured in N) and acts in the direction defined 
by its direction cosines. Multiplying the magnitude of the force by the three direction cosines 
gives the three Cartesian components of the force, F X1 F y and F z , each of which is a scalar, 
and these three numbers (together with the vector’s dimensions) completely specify the vector. 
Multiplying the three components of F by the corresponding unit vectors we can reconstruct 
the original vector by vector addition , 


F = F x x + F y y + F z z, (A.2) 

where the rule for vector addition may be stated as: “Put the 2nd vector so that its tail is at 
the head of the 1st arrow, then the 3rd arrow so its tail is at the head of the 2nd, and the vector 
sum of the three vectors is the vector extending from the tail of the 1st vector to the head of 
the 3rd vector” (see Fig. A.2). Notice that we obtained each of the three vectors that we added 
by multiplying a vector, e.g. x, by a scalar, F x , to get a new vector F x x parallel to the x and 
having magnitude equal to the magnitude of x (1 because it is a unit vector) multiplied by the 
value of the scalar (F x ). 



Figure A.2: The components of F. 
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A.l Dot and cross products 


Multiplying vectors by vectors is a little more complicated as there are different types of mul¬ 
tiplication to consider. We have just now implicitly used the scalar product or “dot product” 
of two vectors, written 

A • B = AB cos 9 (A.3) 

where 9 is the angle between the directions of the two vectors. Taking the dot product of vector 
F with a unit vector gives the component of F in the direction of the unit vector, e.g. 

F x = F • x = F cos 9 X1 F y = F • y = F cos 9 y , F z — F • z = F cos 9 Z . (A.'4) 


The vector product or “cross product” of two vectors, written AxB, is a vector whose magnitude 
is 


|A x B| = ABsinO (A.5) 

where 9 is the angle between the directions of the two vectors, and its direction is given by the 
right hand rule (Fig. A.3). Hence, the Cartesian unit vectors obey the following 

xx y = z, y x z = x, z x x = y, x x z = —y, y x x = —z, zxy = —x, (A.6) 

and xxx = yxy = zxz = 0. We can use these relations to obtain a formula for the cross 
product of two vectors in terms of their components 


A x B = (A x Sc + A y y + A z z) x (B x x. + B y y + B z \ z), 

= 0 + A x ByZ — A x B z y 4 — A y B x z + 0 + AyB z yi + A z B x y — A z B y yi + 0, 


(A.7) 

(A.8) 


x y z 

A X Ay A Z 

B X By B z 


(A.9) 


Two separation vectors which start at the same point, such as A and B in Fig. A.4(a), generate 
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Figure A.3: The cross product and the right hand rule - point your index finger in the direction 
of the first vector (A in this case) and swing your 2nd finger to point in the direction of the 2nd 
vector ( B in this case) and then your thumb is pointing in the direction of their cross product 
(A x B). 


a parallelogram which has an area equal to |A x B|. Three separation vectors which start at 
the same point, such as A, B and C in Fig. A.4(b), generate a parallelepiped. The volume of 
this parallelepiped (also called the scalar triple product ) is 


C • (A x B) = A • (B x C) = B • (C x A). 


(A.10) 


The area of the parallelogram generated by two separation vectors, and the volume of the 
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parallelepiped generated by three separation vectors are useful for constructing the infinitesimal 
surface and volume elements in curvilinear coordinate systems, i.e. Jacobian determinants. 


(a) 



(b) 


height = 

C-n 



Figure A. 4: (a) Cross product of two vectors generates a parallelogram of area |A x B|. 

(b) Scalar triple product of three vectors generates a parallelepiped of volume |C • (A x B)|. 
(n is a unit vector normal to the parallelogram and is in the direction of A x B ) 


A.2 Vector product identities 


a • (b x c) = b • (c x a) = c • (a x b), 
a x (b x c) = b(a • c) — c(a • b). 


B Vector differential calculus 


B.l The Gradient 


(A.ll) 

(A.12) 


In one dimension the differential operator d/dx operating on a function f(x ) gives df /dx, i.e. 
the gradient at x of a plot of / vs. x. In three dimensions the gradient operator V operating on 
a scalar field /(x, y , z) gives a vector field whose magnitude is the steepest gradient, i.e. highest 
rate of increase of / with distance, and is pointing in the direction of the steepest gradient. The 
gradient operator is 


V = 


_ d _ d _ <9 
*Tx +y ¥y +z ¥z 


(B.l) 
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A useful example is the scalar field r n = ( x 2 + y 2 + z 2 ) 71 / 2 , 

Vr n = x^-(x 2 + y 2 + z 2 ) n ^ 2 + y-^-(x 2 + y 2 + z 2 ) n ^ 2 + z-^-(x 2 + y 2 + z 2 ) n ^ 2 , (B.2) 


= x | (x 2 + y 2 + z 2 )^ n/2 l) 2x + y^(x 2 + y 2 + z 2 )( n/2 ^ 2 y 


+ i^(x 2 + y 2 + z 2 )^ 2 - l >2z. 

(B.3) 

n (xx + yy + z z) r n ~ 2 , 

(B.4) 

nrr n ~ 2 , 

(B.5) 

nr n ~ l r. 

(B-6) 


Similarly, for the nth power of the magnitude of the separation vector 
R n = |r — r'| [(x — x') 2 + (y — y') 2 — (z — 2 /) 2 ] 71 / 2 , 

VR n = nR n ~ l R. (B.7) 


We can also differentiate R n with respect to primed coordinates, i.e. using the operator 


_, I ^ d ^ d ^ d 

v = \ x a? + y a^ + z a? 


(B.8) 


with the result 


V'R n = -nR 11 - 1 R. 


(B.9) 


B.2 The Divergence 


The divergence operator is 


V- 


/J ^ d ^d\ 
\*&c +y dy +Z d-z)- 


(B.10) 
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and it operates on a vector field, producing a scalar field which measures the density of the 
whatever is the source of the vector field. The divergence of the vector field A(x, y , z) is 


V-A = 


^ d ^ d ^ i . 
x & +y ^ +z & |A = 


dA x DA,, dA z 


+ 


dx dy 


+ 


dz 


(B.ll) 


As an example, the divergence of A(x, y , z) = 2xx + xyz 2 y + ( x 7 + y + z)z is 


V • A = J“( 2x ) + -§y( X Vz 2 ) + ^(x 7 + y + z) = xz 2 + 3. 


(B.12) 


B.3 The Curl 


The curl operator is 


_ , ^ d _ d ^ d 

Vx = x h _ + y h _ + z h _ x 
ox oy oz' 


(B.13) 
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and it operates on a vector field and produces a vector field which measures the rotation 
(curling) present in the vector field at the point where it is evaluated. In Cartesian coordinates 
its components give the amount of rotation around axes parallel to the x, y and z directions, 


V x A 


x y z 

d_ d_ d_ 
dx dy dz 

Ax Ay A z 


V x A 


dA z dA y \ (dA x dA z \ (dA y 


dA x \ 

dy ) ‘ 


(B. 14) 


(B.15) 


C Vector integral calculus 

C.l Line integrals 


A line integral (or path integral) is an integral of some quantity along some specified path 
in space. A path is one-dimensional and so locations along the path r(A) may be uniquely 
specified by a single parameter, in this case A. As one moves along the path increasing A by dX 
the position vector changes from r to (r + dr) where dr = (dr/dA)dA. For example, the quarter 
circle shown in Fig. C.l (a) may be parameterised in terms of the angle </>, 

r (</>) = acos</>x + a sin</>y. (C.l) 

As one moves along the path increasing (j) by d</> the position vector changes from r to (r + dr) 
where 


dr = 



(^dx ^dy ^dz\ 

i x # +y s +z #j 


dcj) = (—asin^x + acos^y + 0)d</>. 


(C.2) 


The line integral may be of several different types. For example Jp 2 ri |dr| would simply give the 
distance along path V (from ri to r 2 ), and fp 2 ri A • dr would give the integral of the tangential 
component of A along the path T. For example, the line integral of A = frrx around the 
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Figure C.l: (a) A line element dr in a quarter circle in the 1st quadrant of the xy plane, (b) 
Surface element dS in a surface defined by r (u,v) for u\ < u < 112 and v 1 < v < ■ 


quarter circle is 


/ A • dr = / 6(acos0)x • (—asin^x + acos^y) 

/r Jo 


r /2 
— ba 2 / ! 

Jo 


sin 4> cos 4> d(j) = —ba 2 


1 2 

-2 C ° S 


7r/2 

0 


6a2 

T' 


A line integral around a closed loop is represented by jp A • dr. 


(C.3) 

(C.4) 


C.2 Surface Integrals 


The most common type of surface integral encountered in electromagnetism is the flux integral 
which is of the form f s F • dS and gives the flux of the vector field F through the surface. 
For example if F represented the flow velocity (m s -1 ) of an incompressible fluid then the flux 
integral would give the rate at which the fluid pass through the surface (m 3 s -1 ). 

We need two parameters, say u and v to parameterise a surface .S', r = r(it, v). Then the surface 
element can be written dS = A x B, as in Fig. C.l(b), where 


A = 



dr 

B = — dv, 
ov 


(C.5) 
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Hence, for S corresponding to u\ < u < U 2 and v\ < v < v >2 


r ru2 rv2 

/ F • dS = / / F (u, 

JS J u\ J v\ 


dr dr 

—- x — dv du. 
du dv 


(C.6) 


There are two directions perpendicular to any surface, so the direction must be specified - by 
definition it is outward from the volume for integrals over closed surfaces, which are written as 
F • d S. One must always make sure that dS points in the right direction for the problem (i.e. 
into or out of the surface). 

As an example, we shall consider a simple flux integral over the surface of the cone shown in 
Fig. C.2. For points in the surface 

r = rsin# c cos0x + r sin 6 C sin (jry + rcos# c z. (C.7) 

and so it makes sense to use as parameters r (0 < r < R) and cj) (0 < < 27t). 
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Figure C.2: Surface element dS in the surface of a cone. 


Then 


A 


7 —dp = — r sin 9 C sin p dp x + r sin 0 C cos p dp y , 
op 


., <9r 

B = —dr 
or 


sin 6 C cos p dr x + sin 0 C sin p dr y + cos 0 c dr z, 


d S = AxB 


x y 

—r sin # c sin p dp r sin 0 C cos p dp 
sin 6 C cos p dr sin 6 C sin p dr 


z 

0 

cos 9 C dr 


(C.8) 

(C.9) 


(C.10) 


dS — (sin 9 C cos 9 cos p x + sin # c cos 0 sin 0 y — sin 2 9 C z)r dr dp. 


(C.ll) 


We now have the surface element, and for our simple example of a flux integral we will integrate 
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F(r) = F z over this surface, 


r /*27T rR 

F dS= / (Fz) • (sin 9 C cos 6 cos 0 x + sin # c cos 9 sin <J>y — sin 2 9 C z)r dr dq !>, 

J 0 Jo 


2tt /># 


= —F sin 2 9, 


rZ7T r 

Jo Jo 


r dr d<j), 


— —FttR 2 sin 2 9 r . 


(C.12) 

(C.13) 


This is obviously correct since F = Fz is constant, dS -z is everywhere negative and 7r(i?sin0 c ) 2 
is the projected area of the cone perpendicular to z. In this case we didn’t need to parameterise 
the surface or do the integral; we could have just inspected the problem an written down the 
answer! 


C.3 Volume Integrals 


For volume integrals we need three parameters, u,v, and w. The volume element is dV = 
|C • (A x B)| as in Fig. A.4(b) with 


dr dr dr 

C = 7 — du , A = —dv, B = ——dw. 
du dv dw 


(C.14) 


Then 


X Frf3r =/// F \J\ dudv dw 


(C.15) 


where J is the Jacobian determinant, written 


j = dp^z) 

d (l/, 17, w) 


dr f dr dr 
du \dv dw 


(C.16) 


When symmetry permits one can often use cylindrical or spherical polar coordinates for which 
the Jacobian determinants are J — p (cylindrical coordinates p,0,z), an d J — r<2 s i n $ (spherical 
coordinates r,#,</>). 
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C.4 Gradient theorem 


The gradient theorem states 


r r b 


V0 • dr = 0(r 6 ) - 0(r a 


(C.17) 


independent of path, as is easily proved from the definition of the gradient, V</> • dr = dcj). This 
is an exact differential, and so 



4 >(,n) - <Kra)- 


(C.18) 
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C.5 Gauss' theorem 


Gauss’ Theorem (Divergence Theorem) states 


/; 


V • F err = 


j) F-dS 


(C.19) 


where dS is the surface normal which is directed outwards from the volume, and surface S 
bounds volume V. So we could in fact define the divergence of F as 


V ' F = ^o (vf/ dS )’ <C - 20) 

and so the divergence of a vector field F(r) can be thought of as the net outward flux of F per 
unit volume. 


C.6 Stokes' theorem 


Stokes’ theorem (curl theorem) states that 


L 


V x F • dS = 


j) F • dr 


(C.21) 


where loop T bounds surface S', and dr is directed depending on the orientation of dS according 
to the right hand rule for use with Stokes’ theorem as in Fig. C.3(a) — fingers curl around path 
T and thumb points in direction dS. 



Figure C.3: (a) Right hand rule for directions of line integration and surface element when 
using Stokes’ theorem, (b) Three paths joining r\ to r<i- 
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C.6.1 Conservative fields and exactness 

If for a vector field F(r) one finds that V x F = 0, then F is said to be conservative , and Stokes’ 
theorem implies that 


F • dr = 0 


(C.22) 


for any closed loop such as (Ti + 1^) or (Ti + Ts) in Fig. C.3(b). Hence, if F is conservative 


F • dr — (b F • dr — 0 

(Ti+r 2 ) /(ri+r 3 ) 


(C.23) 


Since ri is common to both these closed loops 


r?2 /*r 2 

/ F • dr = / F • dr 

^ri, r 2 


r 2 

ri,r 3 


(C.24) 


and so the integral must be independent of path. This allows us to define a potential 


$(r) = - f F(r') • dr' 
d 1*0 


(C.25) 


such that F(r) = — V4>. Then 


rr rr rr 

- F(r')-dr' = / V$(r')-dr' = / d<$> = $(r)-$(r 0 ), 

Jro J r 0 J r 0 


(C.26) 


and we see that if V x F = 0 then F • dr is exact because V$(r) • dr = d<f> is an exact differential. 


Before moving on, a word of caution is necessary. The result above is true only for simply- 
connected regions, i.e. regions without holes in them. An example of a multiply-connected 
region is shown in Fig. C.4(a). Everywhere along loop Ti we have VxE = 0 and so we might 
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reasonably expect that 


j) E • dr 


0. However, it’s not! Using Stokes’ theorem 


E • dr = f V x E • dS + [ V x E • dS, 

ri Jsi Js 2 

— 0 -\- I V x E • dS, 

Js 2 


= — f E • dr / 0, 

Jr 2 


(C.27) 

(C.28) 

(C.29) 


and this is because the region is multiply-connected and inside we have V x E^O. 
(a) rj (b) rj 




Figure C.4: (a) A multiply-connected region comprising region Si where VxE = 0 surrounding 
a hole (region S 2 ) where VxE/0. (b) As part (a) but with a cut to make region Si simply 
connected. 
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We can make region Si simply-connected by making a cut as shown in Fig. C.4(b). Then, 
integrating around the boundary of the simply-connected region 


® E • dr — 0. 

-'ri+r 3 +r 2 +r 4 

(C.30) 

C.7 Dirac delta function 


Point charges have infinite charge density but finite charge, and we can 
function, S (x — xo), to represent them. The one-dimensional version 6 (x ) 

use the Dirac delta 

is defined such that 

coo 

/ f (x) 6 (x-a) dx = f (a). 

J—oo 

(C.31) 

An important property is 


S(cx) = 

c 

(C.32) 

In three dimensions we have 


<5 3 (r) = 8(x)8(y)5(z), 

(C.33) 

r coo coo coo 

/ J 3 (r)d 3 r = / / / S(x)6(y)6(z)dx dy dz = 1, 

J all space J— oo J— oo J— oo 

(C.34) 

[ f(r)S 3 (r - a)d 3 r = /(a). 

(C.35) 


all space 


In cylindrical and spherical coordinates we need to divide by the appropriate Jacobian 

6 3 (r - r 0 ) = 2 1 S(r - r 0 )8(9 - 0o)5(<j> - <j> 0 ) (spherical), (C.36) 

r z sin a 

6 3 (r - r 0 ) = b(p - po)S((j> - (f>o)8(z - z 0 ) (cylindrical). (C.37) 
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Two important identities are 


V- 



4tt<5 3 (R), 

-4tt6 3 (R), 


where R = (r — r ; ). 


(C.38) 

(C.39) 


C.8 Heaviside step function 


Discontinuous changes of some quantity with position can be represented by the Heaviside step 
function H(x — xq) defined by 


[ 0 (x < x 0 ), 


H(x-xq)=1 1/2 {x = x Q ), 

(C.40) 

[l (x > xq). 


It has the property that 


^ -H(x - x 0 ) = 8(x - x 0 ). 

(C.41) 


An example of its use would be to describe a square voltage pulse of height Vo starting at time 
t\ and ending at £2 for which 

V(t) = V 0 H(t 1 )[l-H(t 2 )}. (C.42) 


Another example would be to describe the electric field on the surface of a spherical conductor 
of radius a carrying charge q. The radial component of the electric field is zero inside and finite 
and radial outside, 

£ - (r) = dh ff<r -“ ) ' <&43) 
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Gauss’ law tell’s us that the charge density is £qV • E 


1 d , 2 x 
= Q~ r ( r E r) > 

(C.44) 

i & ( 2 q tt, C\ 

= E \>Sr\ r 4ueorl H(r “> J ’ 

(C.45) 

= A^dr H{r 

(C.46) 

= 4*r^ (r 0) ' 

(C.47) 


where 5(x) is the Dirac delta function. We know that all the charge resides on the surface 
of a conductor, so that in this case we have a surface charge density a (a) = q/Aiia 2 which 
corresponds to the volume charge density above. 
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D Basics of index notation 


The components of the vector A are denoted where i is its index, i.e. A\ = A x , A 2 = A y 
and A 3 = A z . Similarly, for a matrix, Bij represents the component in zth row and jth column. 


The Kronecker delta is defined 

1 , if i = j 
0 , otherwise 

The Levi-Civita tensor is defined 



(D.l) 


£ijk 


1, if ij k = 123 or 312 or 231 
— 1, if ijk — 1 3 2 or 2 1 3 or 3 2 1 
0 , otherwise 


(D.2) 


Note that is 1 for cyclic permutations of 123, and —1 for anti-cyclic permutations, e.g. 321, 
and 0 otherwise. 


In matrix multiplication [A] = [B] [C] 

3 

Aij = ^ ^ BjfrCkj- (D.3) 

k =1 


In the Einstein summation convention if an index appears twice (dummy index ) then the re¬ 
peated index is summed over for all three values of that index — there is no need to write the 
summation sign — so that 

Aij = BifcCfcj (D.4) 

is entirely equivalent to Eq. D.3. We see that the dummy index k has disappeared while the 
free indices i and j remain. 

In index notation the dot product is 

A B = AiBi. (D.5) 
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Following each summation the number of indices is reduced by one (in the dot product above 
it is reduced from 1 to 0); this is called contraction. Another example of contraction is 

~ Ajk (D.6) 

where this time the number of indices goes from 3 to 2. 

Symmetric tensors obey 

Cij = Cji or C ijk = C k ji (D.7) 

and are unchanged when any two indices are swapped. Anti-symmetric tensors obey 

Dij — Dji or Dijk = Dfcji (D.8) 

and when any two indices are swapped they changes sign. The Levi-Civita Tensor, e is 
anti-symmetric. The Kronecker Delta, 5^-, is symmetric. When you sum over all elements of 
the product of two tensors one of which is symmetric and the other anti symmetric the result 
is zero, e.g. if Aij is symmetric and Bij antisymmetric, then AijBij = 0. 

The gradient operator is 

( d_ d_ d_\ 

~ [dx'dy'dz) ’ 


so that in index notation it is = d/dri. 

The divergence of a vector field D is 

VD = ViDi. (D.9) 

The cross product of F and G, and its zth component are 

F x G = E{jk Gi Fj Gk , [F x G]^ = eijk Fj Gk , (D.10) 

where ei, e 2 and e 3 are the unit vectors in the x, y and z directions. 
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Similarly for the Curl of a vector field H 

V X H — ^ijk ^ jH-ki X H]i — &ijk ^jHk‘ 

A useful identity for the Kroneker delta is 


V iTj 


drj _ x 

dr 3 ~ 


(D.ll) 


(D.12) 


A useful Levi-Civita Identity is 


^ijk^imn 


^jm^kn fijn^km- 


(D.13) 


To make the 1st index of two Levi-Civita tensors the same apply a cyclic permutation to the 
indices of one of them. 
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E Uniqueness of solutions of Poisson's equation 


We will prove by contradiction that the solution of Poisson’s equation is unique. Consider 


the possibility that there could be two different solutions Vi(r) and V^r) for the same charge 
density p(r) and boundary conditions. The corresponding electric fields Ei(r) = — VVi(r) and 
E 2 (r) = -Vb 2 (r) must satisfy Gauss law, 


V • [-Wi(r)] = — and V • [-VF 2 (r)] = ^ 


(E.l) 



V • [V(E 2 - V x )\ = 0. 


(E.2) 


Using product rule Eq. 1.8 we may write 


V • [(E 2 - El) V(E 2 - El)] = V(E 2 - V x ) • V(E 2 - Ei) + (E 2 - Ei)V • [V(E 2 - V)]. 


V • [(E 2 — Ei) V(E 2 — Ei)] = | V(E 2 — Ei)| 2 + 0. 


(E.3) 


Integrating over volume V inside the boundary surface S', 



(E.4) 



(E.5) 


where we have used Gauss’ theorem on the left hand side. 

If V is specified on the boundary surface then on that surface V2(r) = Vi(r), or if E • n is 
specified on the boundary (VV 2 ) • n = (VVi) • n. In either case the surface integral is zero, and 
we are left with 

f |V(E 2 -Ei)| 2 d 3 r = 0. (E.6) 

Jv 

Since the integrand is everywhere greater than or equal to zero, this equation can only be 
satisfied if V(V2 — Ui) = 0, in other words if E 2 (r) = Ei(r). So, we have just have proved that 
the solution for the electric field is unique, and consequently that V^r) = Vi(r) + V cons t where 
Vconst is an arbitrary constant potential. If the potential is specified anywhere on the boundary 
then U 2 (r) = Vf(r) everywhere inside the boundary. 
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F Field line tracing as an initial value problem 


For solving an ordinary differential equation, dy/dx = f(x), given initial conditions y(pc o) = yo, 
one can use Euler’s method with “steps” of length A in x, 

^n+l = %n T ^5 Vn +1 = Vn (F.l) 

Similarly for tracing field lines of vector field F(r) one could use a scheme motivated by Euler’s 
method, as illustrated in Fig. F.l(a) 

r n +i = r n + F(r)A. (F.2) 


A scheme motivated by the 2nd order Runge-Kutta Method for solving initial value problems 
gives greatly improved accuracy, as illustrated in Fig. F.l(b) 


r n+i/2 = r n + F(r n )y; 


r n+1 — r n + F(r n+ 1 / 2 )A. 


(F.3) 




Figure F.l: Comparison of (a) Euler’s method and (b) the 2nd order Runge-Kutta method of 
solving initial value problems as applied to field-line tracing. 

This straightforward technique has been used for producing all the images in this book which 
show electric or magnetic field lines. The best example to show the accuracy of this method is a 
magnetic field line of an infinite straight current, whose field is given by B = (/iqI/ 27vp)<fi. The 
field lines should therefore be circles in the x-y plane, and the two methods discussed above 
are compared for the same step size in Fig. F.2. 
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Figure F.2: Comparison of (a) Euler’s method and (b) the 2nd order Runge-Kutta method 
adapted for tracing the magnetic field of an infinite straight current along the 2 axis. 

G Code for 2D finite-difference method for Laplace's equation 


The IDL code below was used to solve the 2D problem of the potential inside a rectangular 
pipe of square cross section (1 m by 1 m) where the potential is specified on each of the four 
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n_iter=10000 ;;: 

; number of iterations 

nx_cells=100 

ny_cells=100 

nx=nx_cells+2 ;;: 

ny=ny_cells+2 ;;: 

; number of cells in x 
; number of cells in y 
; need extra grid points at cell 
; boundary each side for BC 

x=(findgen(nx)-l)/(nx-3) ;;: 

y=(findgen(ny)-l)/(ny-3) 

; x and y values 

Vx0=0. 

Vxl=0. 

Vy0=0. 

Vyl=l. 

Vsta rt=(Vx0+Vxl+Vy0+Vyl)/4. 
Potential=fltarr(nx,ny)+Vstart 

;;; potentials on the 4 boundary surfaces 

;;; Average of V over boundaries. 

;;; Declare and initialise array for V. 

BoundaryMask=fltarr(nx,ny)+l. 
BoundaryMask(0:1,*)=0. 
BoundaryMask(nx-2:nx-l,*)=0. 
BoundaryMaskC*,0:1)=0. 
BoundaryMaskC*,ny-2:ny-l)=0. 

;;; Set BC mask to 0 at boundary, 1 elsewhere 
;;; After each iteration the potential is 
;;; multiplied by the mask before adding 
;;; the BC potential. 

BoundaryPotential=fltarr(nx,ny) 

;;; Set to BC V at boundary, 0 elsewhere 


BoundaryPotential(0:1,*)=Vx0 
BoundaryPotential(nx-2:nx-1,*)=Vxl 
BoundaryPotential(*,0:l)=Vy0 
BoundaryPotential(2:nx-3,ny-2)=Vyl 
BoundaryPotential(l:nx-2,ny-l)=Vyl 
BoundaryPotential(l,ny-2)=0.5*(Vx0+Vyl) 
BoundaryPotential(nx-2,ny-2)=0.5*(Vxl+Vyl) 


FOR i_iter=l,n_iter DO BEGIN 

;;; Iterate. 

;;; Average over 4 nearby points 


Potential=(SHIFT(Potential,l,0)+SHIFTCPotential,-1,0) $ 

+SFIIFT(Potential ,0, l)+SHIFT(Potential ,0, -l))/4 
Potential=Potential*BoundaryMask ;;; Reset BC 

Potential=Potential+BoundaryPotential 
ENDFOR 

CONTOUR, Potential, X, Y, LEVELS=0. l*findgen(ll) ,/fill 
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sides: U(0, y, z)=V(l, ?/, z)=V(x, 0, z)=0 V, and V(x, l,z)=l V. We choose a grid spacing 
L— 0.01 m, i.e. 100x100 cells (101x101 grid points). After each iteration, the potentials at 
the outer grid points are reset to the potentials on the boundary surfaces. At the corners 
where V changes discontinuously from V—0 to h=l V we set the potential to be the average 
value, i.e. U(0,1, z)=V( 1,1, z)= 0.5 V. The result is contoured in Fig. 2.10. The code uses array 
manipulation methods available in modern programming languages rather than using do loops 
when computing the average potential around each grid point, and in resetting the boundary 
potential. 

The IDL code contains numerous comments, but I will explain some of the procedures called 
which are specific to IDL. 

FINDGEN(n) declares a floating-point array with n elements which have the values: 0,l,2,...n-l. 
FLTARR(n,m) declares a 2D a floating-point array with nxm elements. 

NAME(i:j,*) means the array made up of the elements of the 2D array NAME having the 1st 
index in the rangle i to j and all values of the 2nd index. 

SHIFT (NAME, k,l) means the array constructed from 2D array NAME by shifting the 1st index 
by k and the 2nd by 1. 

CONTOUR, together with some other formatting statements (not shown), contours the result. 


H Dipole term of multipole expansion of vector potential 


The dipole term can be written 



(H.l) 


Using the method of integration by parts, the integrals in the dipole term may be written as 



(H.2) 
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but the 1st integral on the right is zero because d[(r • r , )r / ] is an exact differential. Hence, 


j) (r • r')dr' + 



(r • dr')r', 


(r • dry] . 


(H.3) 

(H.4) 

(H.5) 


Now from the vector triple product rule we have (V x dr') x r = (r-r')dr' 
A.2), and so 

(f (r • r / )dr / = - (j) (r' x dr') x r. 

JVi 2 JTi 


A(r) 


fi o m x r 

47r r 3 ’ 


(r-dr')r' (Appendix 


(H.6) 

(H.7) 


where the magnetic dipole moments for n current loops and a volume current distribution are 



(H-8) 
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I Summary of vector calculus identities 

1.1 Integral theorems 

(i-i) 

( 1 . 2 ) 

(1.3) 

(1.4) 

1.2 Product rules 


Gradient theorem: 


/ V0 • dr — 0(r&) — (j>( r a ) (independent of path). 
J r n 


Gauss’ theorem: 


[ V • F d 3 r = (f F ■ dS (S bounds V). 

Jv Js 


Corollary to Gauss’ Theorem: 


[ V x F d 3 r = - l F x dS (S bounds V). 

Jv Js 


Stokes’ theorem: 


J V x F • dS = j) F - dr (where T bounds S ). 


V(uv) = u(Vv) + v(Vu) (1.5) 

V(A • B) = A x (V x B) + B x (V x A) + (A • V)B + (B • V)A (1.6) 

V • (A x B) = (V x A) • B - (V x B) • A (1.7) 

V • (uA) = (Vu) • A + 7i(V • A) (1.8) 

V x (uA) = (Vu) x A + u( V x A) (1.9) 

V x (A x B) = (V • B)A - (V • A)B + (B • V)A - (A • V)B (1.10) 
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1.3 Second derivatives 


v • (vvo = vV 

(Ill) 

v x (vvO = o 

(1.12) 

V(V • A) = V x (V x A) + V 2 A 

(1.13) 

V • (V x A) = 0 

(1.14) 

V x (V x A) = V(V • A) - V 2 A 

(1.15) 


1.4 Vector operations in Cartesian coordinates (x,y,z) 


dr = dx x + dy y + dz z 
cf' r = dx dy dz 


df ^ df ^ df ^ 


V-A = 


V x A = 


dA x dA n 


dx dy 

dA z dA 
dy 


+ 


dz 


dA 

dz 


X + 


dAx 

dz 


v2 , &L , 

dx 2 <9y 2 dz 2 


9A Z \ _ (dA y dA x \ ^ 

^J y+ {-d^-^i) Z 


(1.16) 

(1.17) 

(1.18) 

(1.19) 

( 1 . 20 ) 
( 1 . 21 ) 
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1.5 Vector operations in spherical polar coordinates (r,0,0) 


dr = dr r + r d6 6 + r sin 9 dcfxj) 
d 3 r = r 2 dr sin 9 d9 d(f> 


df ^ 1 5/~ 1 5/~ 

V/ — "x - r H 7777 0 H : 77 77T </> 
5r r w r sin0 50 


1 5 


1 


9 


V • A = 72 r2A r) + ^( Sin ^) + 


1 dA d 


V x A = 


1 

r sin 0 


5 


r sin0 d9 

dAg 


30 (sin 0.4*)- 


^ l 

r+ - 
r 


r sin 6 d(j) 

1 9A r _ d_ ' 
sin 0 dr T ^ 


9 


1 

+ — 
r 


d dAr 

5r r ^ 50 


0 


1 5 




r 2 dr 


>df 


dr 


1 


d 


r 2 sin 9 d9 


df 


— sin 9-—- + 


d9 


1 g/ 

r 2 sin 2 0 50 2 


( 1 . 22 ) 

(1.23) 

(1.24) 

(1.25) 


(1.26) 

(1.27) 
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1.6 Vector operations in cylindrical coordinates (p,<f>,z) 


dr = dpp + p d(j) <fi + dzz, 


(1.28) 


d 3 r = pdpdcfi dz 


(1.29) 


V/ = 



1 

-77^+7-Z 

p dtp dz 


V A = 


V x A = 


Id, . . 

, 1 

dAj, 

dA z 

~ 7>— \P Ap 

p dp 

)+ ; 

9 + 

d<j) 

dz 

[1 dA z 

&V 


'dA p 

VP <9(/> 

dZ y 

)p+( 

dz 



4> + 



z 


= -—( d/\ Id 2 / d 2 / 
p 9/? \ dp) p 2 dcf) 2 dz 2 


(1.30) 

(1.31) 

(1.32) 

(1.33) 
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J Maxwell's Equations (Static Fields) 


Differential form: 


Vacuum: 


Matter: 



V-E 

P_ 

£ o 

V D 

= PS 

(Gauss’ law) 

V B 

= 0 

V B 

= 0 

(no magnetic charge) 

V x E 

= 0 

V x E 

= 0 

(E is conservative) 

V x B 

= Mo J 

V x H 

= J / 

(Ampere’s law) 


Integral form: 


(j.i) 


Vacuum: 


1 


E • dS = — pd 6 r 
s £ o Jv 


B-dS = 0 


E • dr =0 


Matter: 


® D • dS = pfd 3 r (Gauss’law) 

Js Jv 

B-dS =0 (no magnetic charge) 


E-dr =0 


(E is conservative) 


B • dr = /i 0 / J ■ dS ® H •dr = / Jf-dS (Ampere’s law) 

r Js Jr Js 


(J.2) 

Note: in the 1st row S bounds V; in the 2nd row S is any closed surface; in the 3rd row T is 
any closed loop; in the 4th row T bounds S. 
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